Charged massive particle at rest in the field of a Reissner-Nordstrom black hole 
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The interaction of a Reissner-Nordstrom black hole and a charged massive particle is studied in 
the framework of perturbation theory. The particle backreaction is taken into account, studying 
the effect of general static perturbations of the hole following the approach of Zerilli. The solutions 
of the combined Einstein-Maxwell equations for both perturbed gravitational and electromagnetic 
fields at first order of the perturbation are exactly reconstructed by summing all multipoles, and are 
given explicit closed form expressions. The existence of a singularity-free solution of the Einstein- 
Maxwell system requires that the charge to mass ratios of the black hole and of the particle satisfy 
an equilibrium condition which is in general dependent on the separation between the two bodies. If 
the black hole is undercritically charged (i.e. its charge to mass ratio is less than one), the particle 
must be overcritically charged, in the sense that the particle must have a charge to mass ratio greater 
than one. If the charge to mass ratios of the black hole and of the particle are both equal to one (so 
that they are both critically charged, or "extreme"), the equilibrium can exist for any separation 
distance, and the solution we find coincides with the linearization in the present context of the well 
known Majumdar-Papapetrou solution for two extreme Reissner-Nordstrom black holes. In addition 
to these singularity-free solutions, we also analyze the corresponding solution for the problem of a 
massive particle at rest near a Schwarzschild black hole, exhibiting a strut singularity on the axis 
between the two bodies. The relations between our perturbative solutions and the corresponding 
exact two-body solutions belonging to the Weyl class are also discussed. 
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I. INTRODUCTION 



The problem of the effect of gravity on the electromagnetic field of a charged particle leading to the consideration 
of the Einstein-Maxwell equations has been one of the most extensively treated in the literature, resulting in exact 
solutions (see jlj and references therein) as well as in a variety of approximation methods 0, IMSSSSsS ^] • 

We address in this article the problem of a massive charged particle of mass m and charge q in the field of a 
Reissner-Nordstrom geometry describing a static charged black hole, with mass A4 and charge Q. We solve this 
problem by the first order perturbation approach formulated by Zerilli using the tensor harmonic expansion of 
the Einstein-Maxwell system of equations. Note that a similar approach was used by SibgatuUin and Alekseev 
to study electrovacuum perturbations of a Reissner-Nordstrom black hole. A manifestly gauge invariant (hamilto- 
nian) formalism was developed by Moncrief [11 El [Hi . The relations between these two different treatments of 
perturbations have been extensively investigated by Bicak [T^ . 

Both the mass and charge of the particle are described by delta functions. The source terms of the Einstein 
equations contain the energy-momentum tensor associated with the particle's mass, the electromagnetic energy- 
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momentum tensor associated with the background field as well as additional interaction terms, to first order in m 
and g, proportional to the product of the square of the charge of the background geometry and the particle's mass 
Q^m) and to the product of the charges of both the particle and the black hole qQ)- These terms give origin to 
the so called "electromagnetically induced gravitational perturbation" JT^. On the other hand, the source terms of 
the Maxwell equations contain the electromagnetic current associated with the particle's charge as well as interaction 
terms proportional to the product of the black hole's charge and the particle's mass (~ Qm), giving origin to the 
"gravitationally induced electromagnetic perturbation" . We give in Section II the basic equations to be integrated 
in order to solve this problem. 

Before proceeding to the solution of the most general problem, we recall a variety of particular approaches followed 
in the current literature within the simplified framework of test field approximations 0, 0, hj. In Section III, 
the test field solutions for a charged particle at rest in the field of a Schwarzschild as well as a Reissner-Nordstrom 
black hole are summarized, under the conditions q/m ^ 1, m « and q ^ M, q ^ Q. The case of a test charge 
in the Schwarzschild black hole is reviewed first in Section III A; the corresponding solution is obtained by using the 
vector harmonic expansion of the electromagnetic field in curved space The conditions above imply the solution 
of the Maxwell equations only in a fixed Schwarzschild metric, since the perturbation to the background geometry 
given by the electromagnetic stress-energy tensor is second order in the particle's charge. As a result, no constraint 
on the position of the test particle follows from the Einstein equations and the Bianchi identities: any position of the 
particle is allowed, and the solution for the corresponding electromagnetic field has been widely analyzed ItI, Isj. 

In Section III B we recall this same test field approximation applied to the case of a Reissner-Nordstrom black 
hole. This treatment is due to Leaute and Linet in analogy with the Schwarzschild case, they used the vector 
harmonic expansion of the electromagnetic field holding the background geometry fixed. However, this "test field 
approximation" is not valid in the present context. In fact, in addition to neglecting the effect of the particle mass on 
the background geometry, this treatment also neglects the electromagnetic induced gravitational perturbation terms 
linear in the charge of the particle which would contribute to modifying the metric as well. In Section III C we recall 
how, within this same approximation still neglecting these same feedback terms, Bonnor studied the condition 
for the equilibrium involving the black hole and particle parameters M,q^m as well as their separation distance b. 
He finds the constraint 



If the black hole is "extreme" or "critically charged" satisfying Qj M. — 1, then the particle must also have the same 
ratio qjm — 1, and equilibrium exists independent of the separation. In the general non-extreme case Q/M. < 1 there 
is instead only one position of the particle which corresponds to equilibrium, for given values of the charge-to-mass 
ratios of the bodies. In this case the particle charge-to-mass ratio satisfies q/m > 1. This condition corresponds to 
the overcritically charged case for a Reissner-Nordstrom solution with the same parameters m and q corresponding 
to a naked singularity. 

Already the Leaute and Linet generalization modifies in a significant way the results obtained in the case of a 
Schwarzschild black hole. Indeed a new phenomenon occurs in this case: as the hole becomes extreme an effect 
analogous to the Meissner effect for the electric field arises, with the electric field lines of the test charge being forced 
outside the outer horizon [23| ■ This result by itself justifies the need of addressing this problem in the more general 
Zerilli approach duly taking into account all the first order perturbations. 

The correct way to attack the problem is, thus, to solve the linearized Einstein-Maxwell equations following Zerilli's 
first order tensor harmonic analysis [ij . In Section IIVI we first consider the simplest special case of a massive 
neutral particle at rest near a Schwarzschild black hole. We first show explicitly that a perturbative solution for 
this problem free of singularities cannot exist. We then give the explicit form of the perturbation corresponding to 
a stable configuration when there is the presence of a "strut" between the particle and the black hole, corresponding 
to a conical singularity. This result is a direct consequence of the classical work by Einstein and Rosen j22, who 
considered the Weyl class double Chazy-Curzon solution for two point masses placed on the symmetry axis at 
a fixed distance. We find it very helpful to use a new gauge condition particularly adapted to this problem which 
differs from the corresponding Regge- Wheeler gauge in the Schwarzschild case, by adding additional terms to the 
third Regge- Wheeler gauge condition to take into account the additional terms in the field equations. The resulting 
perturbed metric we obtain is the linearized form of the exact solution representing two coUinear uncharged black 
holes in a static configuration belonging to the Weyl class Q . We review some of their properties and peculiar features 
of these solutions in the Appendix. 

We then turn to the general case of a charged massive particle at rest in a Reissner-Nordstrom background. 
This problem needs the simultaneous analysis of both the vectorial perturbation to describe the electromagnetic 
perturbation and the tensorial perturbation necessary to analyze the perturbation by the particle mass as well as the 
electromagnetically induced gravitational perturbation and the gravitationally induced electromagnetic perturbation. 
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The Einstein-Maxwell equations reduce to a system of 6 coupled ordinary differential equations for 4 unknown func- 
tions determining both the gravitational and electromagnetic perturbed fields. Then, by imposing the compatibility 
of the system (or equivalently, requiring that the Bianchi identities be fulfilled) , we obtain an equilibrium condition for 
the system which coincides with the condition obtained by Bonnor in his simplified approach. This is surprising, 
since our result has been obtained within a more general framework, and both the gravitational and electromagnetic 
fields are different from those used by Bonnor. 

We then succeed in the exact reconstruction of both the perturbed gravitational and electromagnetic fields by 
summing all multipoles. The perturbed metric we derive is spatially conformally flat, free of singularities and thus 
cannot belong to the Weyl class two-body solutions, which are characterized by the occurrence of a conical singularity 
on the axis between the bodies. The asymptotic mass measured at large distances by the Schwarzschild-like behaviour 
of the metric of the whole system consisting of black hole and particle is given by 

A/off = 7W + m + Eint , (2) 

where the interaction energy turns out to be 



-Eint = -m 



b \ b 



-1/2' 



(3) 



Recently Perry and Cooperstock [2^ and Breton, Manko and Sanchez have used multi-soliton techniques to 
study exact electrostatic solutions of the Einstein-Maxwell equations representing the exterior field of two arbitrary 
charged Reissner-Nordstrom bodies in equilibrium. They have shown by numerical methods that gravitational- 
electrostatic balance without intervening tension or strut can occur for non-critically charged spherically symmetric 
bodies in the case of one black hole with Q/Ai < 1 and one naked singularity. However, they have not obtained 
an explicit dependence (in algebraic form) of the balance condition on the separation of the bodies, so we cannot 
compare our analytical formula with their numerical results. Attention is devoted also to the limiting case in which 
Q / M. = q/m = 1. In this case our solution coincides with the linearized form of the exact solution of Majumdar and 
Papapetrou H^H^, where equilibrium exists independent of the separation between the bodies. 

Within our approximation for Q/M. < 1, we find that a horizon seems to form around the particle very close to it 
where the gravitational perturbation is not small compared to the background metric and hence the approximation 
is no longer valid. The study of the full nonlinear problem, namely the corresponding exact two-body solution of 
the Einstein-Maxwell equations is therefore required. Attempts in this direction are being pursued by Belinski and 
Alekseev [27^ by using soliton techniques. 



II. PERTURBATIONS ON A REISSNER-NORDSTROM SPACETIME 



The Reissner-Nordstrom solution is an electrovacuum solution of the Einstein-Maxwell field equations 



Gf_ii, — 87rr™\ 



(4) 



where 



f,P<^F F n F FP'^ 



(5) 



is the electromagnetic energy-momentum tensor. In standard Schwarzschild-like coordinates, the corresponding metric 
is given by 



'f{r)dt^ + fir)-'dr^ + r^{de^ + sin^ ed(j}^) , 



ds' 



(6) 



and the electromagnetic field by 



RN 



-dt A dr 
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where M and Q are respectively the mass and charge of the black hole, in terms of which the horizon radii are 
r± = M ± \/ M?- — Q-^ = ± r. We consider the case IQI < and the region r > r+ outside the outer horizon, 
with an extremely charged hole corresponding to \Q\ — Ai (which implies F = 0) where the two horizons coalesce 
into one. 

Let us consider a perturbation of the Reissner-Nordstrom solution due to some external source described by a matter 
energy-momentum tensor T^i, as well as an electromagnetic current J^. The combined Einstein-Maxwell equations 
are 



(8) 



where the quantities denoted by the tilde refer to the total electromagnetic and gravitational fields, to first order of 
the perturbation: 



rjicra 



FpilFaU ^dfJ.I^FpcyF'' 



Gfiv 



(9) 



note that the covariant derivative operation makes use of the perturbed metric g^i, as well. The corresponding 
quantities without the tilde refer to the background Reissner-Nordstrom geometry and electromagnetic field IZJ. 

Following Zerilli Tl|, we then consider first-order perturbations; the Einstein tensor and the electromagnetic energy- 
momentum tensor @ up to the linear order can be expanded as follows: 



Gni/ — Gill/ ~f" nu 



/T-icm 



(10) 



where 



" 2 ^M^jct 

1 

47r 
1 

An 



Hp-,!') 



F^^Fl", - Ig^^F^^F^^ Kp + ^FP-Fp^h,,, 



(11) 



with fc^ — ^/ict'" and h = ha°'. The expansion of the quantity F^'^.^ appearing in the Maxwell equations turn out to 
be explicitly 



5J^ 



(12) 



where the last term can be interpreted as an effective gravitational current and is given by 

1. 



51^ 



F^'P.^hp" + FP"h%,„ + FPP { kp - -h.,p 



(13) 



As a result, the terms 6x1,^^ - FFh and ST^P 



Ff symbolically represent the electromagnetically induced gravita- 
tional perturbation, while the term <^>/(^) Fh represents the gravitationally induced electromagnetic perturbation 

dill. 

We now specialize this general framework to the case of a perturbing source being represented by a massive charged 
particle at rest. A point charge of mass m and charge q moving along a worldline ^"(t) with 4-velocity U" = dz" /dr 
is described by the matter energy-momentum tensor 



rpflU 

part 



d'''^\x - z{t))UPU'' dT , 



(14) 
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and electromagnetic current 

J^,,, = ^J S^'\x - z{t))U^ dr , (15) 
where the normahzation of the delta functions is defined by 

/.<..(.).'.^i- (16) 

The perturbation equations are obtained from the Einstein-Maxwell equations keeping terms to first order. There 
are two parameters of smallness of the perturbation: the mass m of the particle and its charge q, which can be 
expressed as q = em, where e is the charge to mass ratio of the particle. Therefore, the perturbation will be small 
if m and so q are sufficiently small with respect to the black hole mass and charge, the charge to mass ratio which, 
instead, need not to be small. 

Let us consider the special case of the charged particle at rest at the point r = 6 on the polar axis = 0. The only 
nonvanishing components of the stress-energy tensor (|14|l and of the current density H15|) associated with the charged 
particle are given by 

= £^m'^'Hr~b)s {cose -1) , 

^p°art = {r -b)S [cose -1) , (17) 

since U — f(r)^^/^dt is the corresponding 4-velocity. 

Following Zerilli's [ij procedure, we must now expand the fields h^^^, and f^^, as well as the source terms (tTTI) in tensor 
and scalar harmonics respectively, obtaining a set of first order perturbation equations by linearizing the Einstein- 
Maxwell system . Zerilli showed that in the Einstein-Maxwell system of equations electric gravitational multipoles 
couple only to electric electromagnetic multipoles, and similarly for magnetic multipoles. The axial symmetry of 
the problem about the z axis {0 = 0) allows us to put the azimuthal parameter equal to zero in the expansion, 
leading to a great simplification in the solving procedure. Furthermore, it is sufficient to consider only electric-parity 
perturbations, since there are no magnetic sources 



III. TEST FIELD APPROXIMATION 



Before looking for a solution of the above general problem within the framework of first order perturbation theory 
it is interesting to recall the results in the simplest case of test field approximation, i.e. by neglecting the particle 
backreaction, namely the changes in the background metric due to the mass and charge of the particle. This test field 
approximation requires only the vector harmonic description of the electric field of the test particle alone. 



A. Electric test field solution on a Schwarzschild background 

For any elementary particle the charge to mass ratio is q/m S> 1, so it is natural to address the simplest problem of 
a test charge, neglecting the contribution of the mass, at rest in the field of a Schwarzschild metric, see e.g. the works 
by Hanni Cohen and Wald 0], Hanni and Ruffini '^l, Bicak and Dvorak Linet j^. In standard coordinates 
the Schwarzschild metric is given by 

ds^ = ~U{r)df + Ury^dr^ + r^{de^ + sin^ Odcf)'^) , 

= (18) 
r 

The Einstein-Maxwell equations are 

G^i/ = , 

F^''.^ = 47rJ^ , (19) 

where the only nonvanishing component of the current density associated with the charged particle at rest at the 
point r — h on the polar axis = is given by 

J" = -^(5(r-5)(5(cos6'- 1) . (20) 
27ro^ 
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The electromagnetic stress-energy tensor is second order in the electromagnetic field and can be neglected, so it is 
enough to solve the Maxwell equations in a fixed Schwarzschild background. This solution has been discussed in detail 
by Hanni and Ruffini There is no constraint on the position of the particle in this case. 

The (coordinate) components of the electric field are defined by Ei = FiQ. By introducing the vector potential A^^ 
defined by 

Fa,(3 = 2A^i3,^] , (21) 

which in our case is determined by the electrostatic potential V alone 

Ao = -V , A = 0, (22) 

the electric field becomes 

Er = -Vr , Eg = -Ve , E^ = -V^ . (23) 
The vector harmonic expansion of the relevant components of the electric field is given by 

fIVi 

^. = -E^~^W^"' ^o = -J2^2ir)^ , (24) 



where 



Yio^y^^P,{cos9) (25) 

are normalized spherical harmonics with azimuthal index equal to zero, because of the axial symmetry of the problem 
about the z axis (0 = 0). The expansion of the source term (|20|1 in scalar harmonics is given by 



jO = Y^jO(^r)Y,o, jOir)^^^-^^^S{r-b). (26) 

Therefore, after separating the angular part from the radial one, the Maxwell equations imply 

= £[' + h[ ^^Si + 4nJl> , (27) 
= - fl' , (28) 

where primes denote differentiation with respect to r. By solving Eq. H28|) for £\, and substituting into Eq. (|27|) . we 
obtain the following second order differential equation for E^'- 

= £l" + -8l' - ^$±^£1 + 4vr . (29) 

The knowledge of the function 8\ determines the radial part V/ of the expansion of the electrostatic potential V , from 
Eqs. (123 and HI: 

V ^Y.^i{r)Yia , Vi{r)^£i. (30) 

Putting — fs{ry^^w{r) and making the transformation z = r/M. — 1 Eq. (|29|l becomes 

1 



= (1 - 7J)w" - 2zw' + 



l-z2 



-2V^-^^/2^TT(^^) 5{z-P), (31) 



where primes now denote differentiation with respect to the new variable z and (3 = hj M. — 1. The general solutions 
of the corresponding homogeneous equation are the associated Legendre functions of the first and second kind P^iz) 
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and Q]{z). Using this result and taking into account that Pl{z) = \J — ldPi{z)/dz and Qj{z) = \J z"^ — ldQi{z) / dz, 
Whittaker _,3] and Copson 0| then gave as the two hnearly independent solutions of the homogeneous equation of 129|) 



fiir) 



2'{i + i)inM'+ 

9i{r) = \ 2H\{1 - lyM 
(20! 



(2/ + 1)! _^,_2M)'^Q^^'^'^) 



dr 



(r - 2M) 



dPijzjr)) 
dr 



1 = 0,1,2, 
I = 

; = 1,2, 



(32) 



where Pi and Qi are the two types of Legendre functions. The solution for the electrostatic potential H30|) is then 
given by 



V = qJ2ifi(^)9i{rWb - r)+gi{b)fi{r)d{r - 6)]P,(cos0) . 
I 

Copson Q showed that this solution can be cast in closed form using the following representation formula 

xt ~ cos Q 



(33) 



[x2+t2-2a;<cos6i-sin2 6I]V2 



x^(t - x) + W{x -t)- {x^ - l)(i2 - 1) 



E 

1=1 

dPi{x 



21 + 1 



1(1 + 1) 



dQi{x) 



dx 



dx 

dQijx) 
dx 



dPijx) 
, dx 



■dit - x) 



i){x - t) 



Pi (cos 0) 



(34) 



However, the solution he found did not satisfy the boundary conditions; the corrected version was presented about 
fifty years later by Linet 0: 



Vs = 



q [r- M){h- M)- W^cosO qM 



br 



br 



with 



Ds = [{r ^ M)^ + (b- M)^ -2{r - M)ib- M)cos9 - M^sin^ 9]^/^ 
The components of the electric field are then easily evaluated 



(35) 
(36) 



Er 



q 

Jjj.2 



M 



M{b -M)+M'^ cos d r[{r - M)ib - M) - cose][{r - M) - (b ^ M) cost 



Ds 



Ee = qbfMrUr) 



sin 6 



(37) 



The properties of the above solution has been analyzed in detail by Hanni and Ruffini ■ They derived the lines of 
force by defining the lines of constant flux, and also introduced the concept of the induced charge on the surface of the 
black hole horizon, which indeed appears to have some of the properties of a perfectly conducting sphere terminating 
the electric field lines. 



B. Electric test field solution on a Reissner- Nordstrom background 

The generalization of the above treatment to the case of a charged test particle at rest near a Reissner-Nordstrom 
black hole was discussed by Leaute and Linet the Schwarzschild case, the problem is solved once the 

solution of Eq. (|29|1 with fs{r) < — > f{r) is obtained. The corresponding solution for the electrostatic potential is 
again of the form H33|) with functions 

, _ (2^ + 1)! {r-r+){r~r-)dQi{z{r)) , ni9 

' " 2'(Z + l)!l!r'+i r dr i-U,L,z,... 

1 1 = 

5/(0 = \ 2H\{1 - ly.r' {r - r+){r - r_) dPi{z{r)) (38) 

i — 1, 2, ... 



(20 ! r dr 
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where now z = (r — M)/T. Leaute and Linet showed that this solution can be cast in the closed form expression 

q {r~M){b-M)-T^CO^e qM 

^^^ = 6;: D^^ + 17' 

with 

i^RN = [{r-Mf ^{b-Mf -2{r-M){h - M) cose- T'^sm^eY''^ . (40) 
The components of the electric field of the test particle alone are then easily evaluated 

^ J A. M{h- M) + V'^cose r[{r - M){h - M) -V^ cos e][{r-M)-{b~M) cost 

hir — -nr { M. ' 



br^ I i?RN 
sin 6 



Eg = qbf{b)rf{r)-^ . (41) 
The black hole has its own electric field and electrostatic potential 

i?f H = ^ , 1^^" = ^ . (42) 

This test field approach has been largely used in the current literature |0, 0, . It is interesting to analyze 

what conceptual differences are introduced in the properties of the electric field of the test particle in the context of 
a charged Reissner-Nordstrom geometry compared and contrasted with the Schwarzschild case. The concept of the 
induced charge on the horizon and the behavior of the electric lines of force have been addressed in . We recall here 
the main result: as the hole becomes extreme an effect analogous to the Meissner effect for the electric field arises, 
with the electric field lines of the test charge being forced outside the outer horizon. We are witnessing a transition 
from the infinite conductivity of the horizon in the limiting uncharged Schwarzschild case to the zero conductivity of 
the outer horizon in the extremely charged Reissner-Nordstrom case. 



C. Bonnor's equilibrium condition 

In this same test field approximation Bonnor |19| | has addressed the issue of the equilibrium of a test particle of 
mass m and charge q at rest at r = &, = outside the horizon of a Reissner-Nordstrom black hole. By considering 
the classical expression for the equation of motion of the particle 

mC/"V„f7'3 = qPf^^U^" , (43) 
with 4-velocity U" — /(»')^^^^^o i found the following equilibrium condition for such a system 

From this equation it follows that there exist equilibrium positions which are separation-dependent, and require either 
< m} and > A^^ or > vr? and < . A condition which is sufficient, but not necessary for the equilibrium 
is |g| — m and \Q\ — M.\ it represents a special case of the Newtonian condition qQ = mM.^ so the equilibrium can 
occur at arbitrary separations. We will again address this issue on the existence of equilibrium conditions in Section 
Ivl using the correct treatment which takes into proper account the perturbation of the Reissner-Nordstrom field 
induced by the charge and mass of the particle. 



IV. PERTURBATION ANALYSIS: THE SCHWARZSCHILD CASE 



Let us consider now the problem consisting of a neutral particle of mass m at rest on the polar axis near a 
Schwarzschild black hole in the framework of first order perturbation theory. The Einstein equations are 

G^u = 8^r;r , (45) 

where the perturbed Einstein tensor is defined in ©, and the source term is given by ()17|l . We need to expand the 
Einstein tensor as well as the source term in tensor and scalar harmonics respectively, obtaining a set of first order 
perturbation equations, once a gauge is specified. 
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The static perturbations due to a point mass in a Schwarzschild black hole background has been studied by Zerilli 'sj 
in the dynamical case following the Regge- Wheeler [s^l treatment. The geometrical perturbations /i^^ corresponding 
to the electric multipoles are given by 

dYio 



sym e-'''H2Yi^ 



ho 
hi 



86 
dYio 
86 



sym 



sym 



sym 



sym 



KYio + G 



sym 



d'Y, 



10 



86^ 







r2 sin^ei ( KYiQ + G cot 



,8Y, 



10 



89 



(46) 



where the symbol "sym" indicates that the missing components of hf^^, are to be found from the symmetry h^u — h^^, 
and e"^" — fs{r) is Zerilli's notation. 

This general form of the perturbation can be simplified by performing a suitable gauge choice. Consider the 
infinitesimal coordinate transformation 



(47) 



where the infinitesimal displacement ^'^ is a function of and transforms like a vector. This transformation induces 
the following transformation of the perturbation tensor hf^i^: 



7 (now) „ J -2f/ ^ 



(48) 

where the multipole expansion of the second term can be easily obtained by expanding in the electric-type vector 
harmonics 



8Yi 

AoYio8t + AiYiodr + A2^de 

86 



(49) 



The new metric perturbation functions are then given by 



H, 



(new) 




Ho 



2M 



r{r - 2M) 

2M 
r{r - 2M) 



Ai - 2Ai' , h^r^^ = hi- 



1+1 Mo 



r J 
r 

r -2M 



(new) 
^0 



Ai - r^Aa' , 



/lo + ( 1 - ^ ) , 



^(ncw) ^K--Ai , G("°") G - 2^2 



(50) 



where a prime denotes differentiation with respect to r. 



A. The Regge- Wheeler approach 

We use the Regge- Wheeler [s^ gauge to set 

^(RW) ^ ^(RW) ^ ^(RW) ^ 

This specialization is accomplished through the gauge functions 



Ao^-^^ho, A,= [l-—)[^ki--G') , M^-, (51) 

so that the metric perturbation functions (|50|l expressed in the Regge- Wheeler gauge as combinations of metric 
perturbations expressed in an arbitrary gauge are given by 



H r{r - 2M) 



^H, + r{r - 2M)G" + (2r - 3X)G' - 2 f 1 - ^ ^i' - ' 



^(RW) =K+{r- 2M) ( G' - ^/u ) . (52) 
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The general perturbation (|46|l thus becomes 



IIVI 





HiYio e-''H2Yio 
r^KYio 



















r2 sin^ OKY,. 



10 



(53) 



where superscripts indicating the chosen gauge have been dropped for simpHcity. 

Therefore, the independent first order perturbations of the quantities appearing in the Einstein field equations 145|) 
are given by 



Gqo 
Gil 

G12 

Gqi 



1 
2 

26"" 



2K" - -H2' 
r 



6 



K' -2 



1 

- + — 



(i^o + H2) 



iX + l)H2-Ho + XK] Wio 



1(2 , 



+ -{H^-H2 



^^[ff2 - (A + l)Fo + \K] \Ym , 



(54) 
(55) 



V,' + -]K' - Ho' 



^ + -)^2'-(^ + -)Ho' + 2(A + l)V(^o-i?2)h 



d'Yio 



Ho' + K'-[^ + -]H2 



'-}H 



\dYio 



A e" 
— + - 



G, 



02 



^part 



r 

— {Hi +'^sHi} — 



Hi ^Yia 



16n 



-Af^nY 



00 ^iO 7 



where X — ^ {I — 1) {I + 2) and the quantities with the subscript/superscript s stand for 



fs(r) = 1 



2M 



A' 

JT-i 



00 



■my/21 + 1 



(56) 

(57) 

(58) 

(59) 
(60) 

(61) 



The angular factors containing derivatives vanish for / = 0; moreover, the two angular factors in the expression l|89|) 
for G22 are not independent when Z = 1 (in fact, d^Yio/dO^ ~ — ^lo)- Therefore, the cases ^ = 0, 1 must be treated 
separately. 

For all higher values of /, the Einstein field equations 145|l imply that the corresponding curly bracketed factors on 
the left and right hand sides are equal, so that the system of radial equations we have to solve is the following: 



2K" - - W' 
r 



K' -A 



1 vJ 



W 



{W + K)+Alo , 



'-W' 



2\p^''^ 



^ K" 



K' - W" - 2\ v^' + - ]W' , 



-W + K' - Vs'W 



(62) 

(63) 

(64) 
(65) 



since 



Ho = H2 = W , Hi = 



(66) 



We are dealing with a system of 4 ordinary differential equations for 2 unknown functions: K and W. Compatibility 
of the system requires that these equations not be independent. The second order equations (|62|l and H64|) show that 
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the first derivatives of both the functions K and W are not continuous at the particle position with the same jump 

[K'] = [W] = -A^^V^TTifM-^'^ , (67) 

where the quantity [K'] = K'j^ — K'_, K'^ — K'{bzt e) is obtained by integrating both sides of Eq. Ht)2l) with respect 

to r between b — e and b + e and then taking the limit e ^ 0. Only the pair of first order equations and (jH5|l thus 
remain to be considered. 

It is useful to introduce the the new combinations 

X = K -W , Y = K + W , (68) 

so that the system we have to solve for the unknown functions X and Y is the following 

= --X' - ^(X' + Y') + -^^X , (69) 

= X'- ^(X-r) . (70) 

Solving Eq. (|7U|I for Y and then substituting it into Eq. (|69|l gives the following equation for X 

= r(r - 2M)X" + 4(r - M)X' - {I + 2){l - 1)X . (71) 
Putting X ~ fs{r)^^/^w{r)/r and making the transformation z = r j M. — 1 the previous equation becomes 

1 



= (1 - z^)w" - 2zw' + 



w , (72) 



where primes now denote differentiation with respect to the new variable z. The general solutions are the associated 
Legendre functions of the first and second kind Pl{z^ and Q\(z\ implying that 

X = (r- 2X)-i[ci/K0 + C2gi(r)\ = ci^i + c^X^ , (73) 

where the functions fi{r) and gi{r) have been defined in H32|) . Choosing the arbitrary constants ci and C2 in order 
that the function X be continuous at the particle position and to satisfy regularity conditions on the horizon and at 
infinity implies that the solution can be written as 

X =Afx[Xi{r)X2{b)d{b-r)+Xi{b)X2{r)^{r~b)] , (74) 

where J\fx is an arbitrary constant. The corresponding solution for Y can be easily obtained from Eq. H7(J|I : 

Y = X + -^{r-2M)X' . (75) 
M 

Inverting the relations H68|l yields immediately the solutions for the functions K and W: 

K^^ + ^i^-^^)^' ' W = -^{t~2M)X' . (76) 

The value of the arbitrary constant Mx is determined by imposing the condition H67|l : 

Mx = 8V^^A^V2]TT/s(6)-'/'[X2"(6) - Xi"(fe)]-i . (77) 

Next consider the case / = 0. The relevant equations come from quantities (|54|I - H60|) which do not contain angular 
derivatives: 

= r(r - 2M)K" + (3r - hM)K' - (r - 2M)H^ + K - H2 + - (78) 

= (r - 2M)Ho' -{r- M)K' + H2 - K , (79) 
= r{r - 2M)[K" - Ho"] + {r - M)[2K' - H2'] - (r + M)Ho' . (80) 
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It is easy to show that looking for solutions of the form Hq = H2 = W leads to a system of equations which coincide 
with Eqs. I|t)2|) ~ (|64|l for / = (or A = —1). An analogous situation occurs in the remaining case I — 1: the relevant 
equations coming from quantities 154|l -16() () 



= r(r - 2M)K" + (3r - bM)K' - (r - 2M)H2' - 2H2 + -- — ^777^^0 : (81) 



1 

2r-2A^' 

= (r - 2M)Ho' - (r - M)K' + H2 - Hq , (82) 
= r(r - 2M)[K" - W] + {r - M)[2K' - H2'] - (r + M)Ho' + Hq - H2 , (83) 
= r(r - 2M)[K' - Hq'] - (r - M)H2 + (r - ?,M)Hq , (84) 

reduce to the system (|62|l " (|65|) where we set Z = 1 (or A = 0). Therefore we can take the solutions H76() for all values 
of Z. 

The reconstruction of the solution summing over all multipoles is nontrivial. Furthermore, the first derivatives 
of both the functions K and W are not continuous at the particle position implying that the perturbed Riemann 
tensor is singular there. Indeed, a singularity-free solution for this problem is obviously impossible, since there is no 
external force to oppose the infall of the particle towards the black hole, so that equilibrium cannot be reached in 
any way. We will see in Section IV B how adopting a gauge different from the Regge- Wheeler one gives rise to a 
more convenient form of the gravitational perturbation functions, yielding a closed form expression for the perturbed 
metric by summing over all multipoles. The singular character of the solution will be manifest in this new gauge. 

Of course, one can get regular solutions even in the Regge- Wheeler gauge either by modifying the symmetry of the 
problem e.g. by the introduction of angular momentum or transverse stresses to balance the gravitational attraction 
making the resulting configuration as stable. Consider, for instance, the solution describing a thin spherical shell of 
matter with suitable isotropic pressure at rest around a (concentric) Schwarzschild black hole. An exact solution for 
this problem has been found by Frauendiener, Hoenselaers and Konrad |33| | . The spherical symmetry of the problem 
(with only two additional components Teg and T^^ of the stress-energy tensor of the source, while Trr is assumed 
identically zero) allows to get a regular solution, to which only the monopole / = contributes. A problem with less 
symmetry would require the addition of anisotropic transverse stresses, and all higher values of I may contribute in 
this case. 



B. The perturbation analysis of the Weyl class 

As it is well known, there exist exact solutions to the vacuum Einstein's field equations representing the nonlinear 
superposition of individual static gravitating bodies in an axially symmetric configuration: the solutions belonging 
to the Weyl class (see e.g. ,1,]). These solutions are not singularity- free, but exhibit singular structures as "struts" 
and "membranes" necessary to balance the bodies. We refer to the Appendix for a summary on these many-body 
solutions belonging to the Weyl class. In particular, we are interested in the exact solutions corresponding to a point 
particle at rest above the horizon of a Schwarzschild black hole as well as to a pair of coUinear Schwarzschild black 
holes. We thus expect to find a solution of the first order perturbation equations obtained following Zerilli's approach 
that represents just a linearization of these known exact solutions. 

Let us start again with the general form of the perturbed metric (|46f) : first of all, we use the available gauge freedom 
to eliminate two of the off-diagonal terms, corresponding to the first two Regge- Wheeler conditions 



= hi — Q 



(85) 



We then have 



IVII 



HiYio e-^^H2Yio 







2 i KYio + G 








r2sin2 0(xy,o + Gcot-^^'° 



961 



(86) 



The remaining gauge freedom will be exploited later. 
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The independent first order perturbations of the quantities appearing in the Einstein field equations H45|) are 
Goo = — 



2K" ~ 2(A + 1)G" - -H2' + (W + -][K'-iX + 1)G'] 
r V r ,' 



1 l^s 

1^ — 

,2 



r 



{Ho + H2 



26"= 



[(A + 1)H2 -Ho + \K] Wio , 



(87) 



Gil 
G22 

G12 
Goi 
G02 



rppart 
-'OO 



2 r 



2e- 



' + - Wf'-ifo" - + - - P^ + - Uo' 1/0 



[i/2-(A + l)i/o + Ai^]^r/o 



iv' 1 



cot 6* 



H2 - r^e"' 


G" + 















BY, 



10 



80 



2 r ' " 



1^ 
J 09 



A 

— + - 



1 

f - 

r 

dYio 



Hi }Yio , 



(89) 

(90) 

(91) 

(92) 
(93) 



where the gravitational source term coming from the point particle is still given by Eq. (|61() . The angular factors 
containing derivatives vanish for I = 0; moreover, the two angular factors in the expression (|89|l for G22 are not 
independent when / = 1. Therefore, the cases / = 0, 1 will be treated separately. 

For all higher values of /, the Einstein field equations H45|l imply that the corresponding curly bracketed factors on 
the left and right hand sides are equal, so that the system of radial equations to be solved is the following: 







„2l^B 



6 



2K" - 2(A + 1)G" - -H2' +{j^s' + -][K' -(X + 1)G'] - 2 { — + ^ ] {Ho + H 



1 I^s 



-^[(^ + 1)^2 -Ho + XK] + Ala , 
- lHo'-(iys' + -) [K'-{X + 1)G'] + 



= K" 



2e- 



-[H2-{X+l)Ho + XK] 



= H0-H2- r^e 



-Ho' + K'-G' 



l]K'-Ho"-{^-^ + -]H2'-{^ + -]Ho' 

2„i's 



G" + 













(94) 
(95) 
(96) 
(97) 



2 r 



Ho 



Ho 



since Hi = 0. 



1. The new gauge condition 



Consider first the case I > 2. We are dealing with a system of 5 ordinary differential equations for 4 unknown 
functions: Ho, H2, K and G. Compatibility of the system requires that these equations not be independent. It is 
easy to show that Eq. H96|l is identically satisfied by substituting into it the quantities K' , K" and H2 obtained from 
Eq. H98|l . the same equation after differentiation with respect to r and Eq. (j^Zjl, respectively. 

By imposing the following further gauge condition on the perturbation functions 

Ho + H2 = 2[K-{X + l)G] , (99) 

instead of the third Regge- Wheeler gauge condition G = 0, we show in the following that with this specific choice 
we can obtain in a manifest way the linearized form H146|) of the metric of two collinear Schwarzschild black holes 
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belonging to the Weyl class. In fact, direct inspection of the metric components (|147|) below clearly shows that they 
certainly satisfy the above condition, once expanded in multipoles according to the general form H46|l. It is clear that 
this remarkable result cannot be obtained by using the Regge- Wheeler gauge, since in that case the perturbation 
analysis remains within the framework of a single Schwarzschild solution. We will refer to this new gauge as BGR 
gauge. 

Summarizing, our gauge choice is completely specified by the three conditions 

^(BGR) ^ ^(BGR) ^ Q ^(BGR) ^ ^(BGR) ^ 2j^(BGR) _ ^ I)g(BGR)] . (IQO) 

The metric perturbation functions (|50|l expressed in this gauge are thus given by 

rriBGR) _ rr , A fr(BGR) _ rr , u h ' 

ijf ^ 2(A + 1)(2A, -G)-Ho + 2K- , 

r r — 2M 

^(BGR) ^K--A^, G'^GR) = G - 2^2 , (101) 

r 

since 

Ao = -y^h, , (102) 
while the remaining gauge functions must satisfy the following equations: 

Ai'^-^—^A,-{\ + l){2A,-G) + ^{Ho + H2)-K, A^' = - \ +^ . (103) 
rr — 2M. 2 r{r — 2M) r'' 

Differentiating both sides of the second equation above with respect to r and using the first equation in it leads to 

^2" = -r^^T7T^2+A, A=f^-2^-—^-—[2{X+l)G + Ho + H2-2K], (104) 
r{r-2M) H 2r(r-2A^) 

which involves the function A2 only. 



2. Explicit solution for the Weyl-type perturbations 

Let us look for a solution of the system consisting of equations , , and , taking into account the 
condition (|99|) . By solving the gauge condition (|99l) for G, and substituting it into all the other equations, we have: 

r(r - 2Mf[H2" + Hq"] + (r^ - 'iMr + 2M^)H2' + (Sr^ - llXr + \QM^)H^' 

-{r - 2M)[{f + I + 2)H2 -{1 + 2){l - l)K] - Al^ , (105) 

= (r - M)H2' -{r- 3M)Ho' - 2H2 + l{l + l)Ho - {I + 2){l - 1)K , (106) 

= (r - 2M)[2K" - H2" - Ha"] + (r - M)[2K' - H2' - + H2 - Ho , (107) 

= r(r - 2M)[H2' -{f + l- l)i?o' + {I + 2){l - l)K'] - l{l + l)[{r - M)H2 - {r - 3M)Ho] . (108) 

Next we eliminate the function K by solving Eq. (|108|l for K', substituting it and its derivative K" into Eq. (|107|l . 
and then by solving Eq. H106|l for substituting it into Eq. (|105|l . thus obtaining a pair of equations involving the 
functions Hq and H2 only 

= r(r - 2M)[Ho" + H2"] + 4(r - 2M)Ho' - l{l + l)[iJo + H2] + ^^^^00 , (109) 
= r(r - 2M)[H2" - Hq"] - iMHo' - l{l + l)[H2 - Hq] . (110) 
By subtracting the previous equations, we obtain the following second order differential equation for the function Hq'- 

= r(r - 2M)Ho" + 2(r - M)Ho' - l{l + l)i7o + 2{/-2M) ^^^° ' ^^^^^ 
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Making the transformation z = r / M — \ leads to 



= (1 - z'^)Ha" - 2zHo' + l{l + l)-ffo - V2r+T 



/3-1 
13+1 



1/2 



5{z - P) , 



(112) 



where primes now denote differentiation with respect to the new variable z and (3 ~ b/A4 — 1. The general solutions of 
the corresponding homogeneous equation arc the two types of Legendre functions Pi and Qi. After imposing regularity 
conditions on the horizon and at infinity, the solution of Eq. (|112|l is then given by 



Ho - 4V^V2rTT— [Pi{z)QiiP)^{b -r)+ F,(/3)Q,(z)^(r - b)] . 



(113) 



From Eq. H11U|) . the corresponding solution for H2 is given by 



MJ 

dQi{z{r)) 



{b-M)fs{hY'\{r-'^M) 



dr 



r—b 



dr 



dr 



dQijzjr)) 
dr 



d{r ~ b) 



(114) 



Expressions for the remaining functions G and K can be easily obtained from Eq. H97|) and the relation (|99|l respec- 
tively. 

Finally we have to analyze the cases 1 — 0,1 separately. Consider first the I = case. The relevant equations come 
from quantities (|87|I - H9UI) which do not contain angular derivatives. They now have the form 



= r(r - 2M)K" + (3r - 5M)K' - (r ~ 2M)H2' + K ~ H2 + 2(^r ~ 2M) '^°° ' 
= (r - M)K' - (r - 2M)Ho' + K - H2 , 

= r{r-2M)[K" ^Ho"] + {r-M)[2K' -H2']-{r + M)Ho' . 
The algebraic gauge condition then becomes simply 

Ho + H2 = 2K . 



(115) 

(116) 
(117) 



(118) 



Thus the function G drops out of the equations and can be chosen to be equal to the corresponding one for I > 2 
evaluated at I = without loss of generality: in fact inspection of the perturbed metric tensor (|151|l shows that the 
corresponding angular factors vanish identically for ^ = 0. 

Then by solving Eq. (|118|l for K, and substituting it into equations (|115|I - H117|) . we obtain the following three 
equations involving the functions Hq and H2 only: 

= r(r - 2M)[Ho" + H2"] + (3r - 5M)Ha' + (r - M)H2' + Ho - H2 + ^r^^^oo . (US) 

= (r - M)H2' - (r - 3M)Ho' + Ho - H2 , (120) 
= r(r - 2M)[H2" - Ho"] - AMHo' . (121) 

By subtracting the last two of the previous equations from the first one, we obtain the following second order differential 
equation for the function Ho'. 



= r(r - 2M)Ho" + 2{r - M)Ho' + 



which coincides with the Hlll|l for I = 0. So, its solution is simply 



2M 



2{r-2M) 



-A 



00 ' 



Ho = -2V^^Mby/' 
from Eq. I|117(l . H2 results to be given by 



In 1 — i9(6-r)+ln 1 h9(r - 6) 



2M 



H2 - (3-24-)i/o + 8V^^/.(6)i/2 



MJ 



{r - M){b - M) - 
b{b-2M) 



■d{b~r)+-d{r -b) 



(122) 



(123) 



(124) 
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while the function K can be easily obtained from relation H118|l . 

Only the case / = 1 remains to be considered. The two angular factors in the expression H89|l for the G22 component 
of the Einstein tensor are not independent when Z = 1, since Yio = cos0, so the two curly bracketed factors must be 
considered together. The relevant equations coming from quantities (|87|) " (|90|l are thus given by 

3 

= r(r - 2M)[G" - K"] + (3r - bM)[G' - A^] + (r - 2M)H2' + 2H2 - 2(7^2A4)^°° ' ^^^^^ 

= (r - M)[K' - C] - (r - 2M)Ho' + Ho - H2 , (126) 
= r(r - 2M)[K" - G" - Ho"] + (r - M)[2{K' - G') - H2'] - (r + M)Ho' + Hq - H2 , (127) 
= r(r - 2M)[K' - G' - Ho'] - (r - M)H2 - (r - iM)Ho . (128) 

The gauge condition (pl^ becomes 

Ho + H2 = 2{K - G) . (129) 

By solving Eq. p29|l for G, and substituting it into Eqs. (|125|I - H128|I . we obtain the following fom- equations involving 
only the functions Ho and H2'. 

3 

= r(r - 2M)[Ho" + H2"] + (3r - bM)Ho' + (r - M)H2' - AH2 + y^^M^lo , 

= (r - M)H2' ~{r- 3M)Ho' + 2{Ho - H2) , 

= r(r - 2M)[H2" - Ho"] - AMHo' + 2{Ho - H2) , 

= r(r - 2M)[Ho' - H2'] + 2{r - M)H2 - 2(r - 3Al)i/o • (130) 

By subtracting the second and third of the previous equations from the first one, we obtain the following second order 
differential equation for the function Ho 

3 

= r(r - 2M)Ho" + 2{r - M)Ho' - 2Ho + ^(^^^Al) ^"" ' ^^^^^ 
which coincides with Eq. (|lll|l for I — 1] therefore its solution is simply 



Ho = 4x/3^—/, (6)1/2 [p^(^)Q^(^)^(& _ + p^{^i3)Q^(z)§{r - b)] , (132) 



where 



Pi{z) = z, Qi(^)^f ln(^i±I^ -1 . (133) 

Furthermore, the third of Eqs. (|130|l coincides with Eq. H11U|) for ^ = 1; the corresponding solution for H2 is then 
given by Eq. I|114|l evaluated at / = 1. The function K remains undetermined, and can be chosen to be equal to the 
corresponding one for I > 2 evaluated at ^ = 1 without loss of generality. Then the function G can be easily obtained 
from Eq. ifT^ . 

3. The analytic solution summed over all values of I 

We want now to reconstruct the solution for the gravitational perturbation functions Hq, H2, K and G for all values 
of I. We will denote by a bar the corresponding quantities summed over all multipoles. 
Consider first the solution (|113|l for Ho- The sum over all multipoles turns out to be 



Hq — ^ HoYio 



2^/s(&)'/' HoYio [Pi{z)Qiip)il}{b -r)+ Pi{fi)Qi{z)d[r - b)] P,(cos( 



^^JM" , (134) 



(=0 
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where the quantity Ds is defined in (|36|l and the following representation formula has been used: 



Y,{'^l + ^)[Pi{x)Qi{tW-x)+Pi{t)Qi{x)^{x-t)]Pi{cose) . (135) 



[a;2+t2-2a;icos6'-sin2 6']i/2 
Consider then the solution pi4|l for H2. The sum over all multipoles turns out to be 



oo oo ^ 

1=0 1 = 1 



2 /=0 



00 



!=0 



E 



2/ + 1 



dQi{z{r)) 



dr 



dr 



[r - M){b - M) - M'^ 



b{b-2M) 



- 2— f 

" Ds^'^' b{b^2M) 



fs{b) 



1/2 



1 - 



dr 

d{b - r) + d{r - b) 

r - M — {b — M) cos ( 
D~s 



dQijzjr)) 
dr 



?9(r - b) 



Pi (cos ( 



(136) 



where z = r/A^ — 1 and the representation formula H34I) has been used. 

In order to find the sum over all multipoles of the remaining gravitational perturbation functions K and G it proves 
more convenient to proceed as follows, rather than deriving first the corresponding multipolar solutions and then 
trying to sum the series. Consider Eqs. (|97|l and (|98|) : since the parameter I does not appear explicitly, they remain 
valid for the corresponding summed quantities as well and can be rewritten as 



= Ho-H2-dr [r{r - 2M)drG] , 

= r{r-2M)dr[-Ha + K -G]-{r-M)H2 + {r-3M)Ho 



(137) 
(138) 



These equations must now be treated as partial differential equations rather than ordinary differential equations. In 
fact summing over the spherical harmonics leads to the barred functions which are thus depending also on the angular 
variable 9. After integration on the radial variable each function G and K will be determined up to an arbitrary 
function of 9, which will be then chosen in order the perturbed metric to satisfy the Einstein field equations H45|l. 
Substituting the solutions H134|) and (|136l) into Eq. 1137|) gives the following first integral 



drG 



iMm ^ ^^y/2 r - Ds + ki{9) 



b{b-2M)-' ' ' r[r-2M) ' 
where k\ (ff) is an arbitrary function of the polar angle. A further integration gives 



(139) 



G 



b{b-2M) 
cos 9 In 



1/2 



-4(/32 



-In 



1) sin^ 



(3 cos ( 



J 



-/31n 



1 



2J{z(3 -cos9 + J) 



2J{P~ ZC0S9 + J) 
z2 - 1 



(/32 - l)(z2 - 1) 

fci(6i)arctanh(z) + ^2(6*) 



(140) 



where z — r/Ai — 1, f3 — b/Ai — 1, J — Ds/M and k2{9) is another arbitrary function of 9. The solution for K 
follows immediately by integrating Eq. H138|) : 



K = Ho + G 



AMm 
b{b^2M) 



fs{by/'\-ln 



z — (3 cos 6 + J 



Vi2^ 



-f3\n 

2^ 



2 J{zl3 -cos 9 + , J) 



(/32 - l)(z2 - 1) 



+fc3(^) , (141) 



where ^3(0) is a third arbitrary function. 

Therefore, starting from the expansion (|46|l for the gravitational field, the perturbed metric summed over all 
multipoles turns out to be 



dS^ = -e^^ll - Ho]dt^ + €-"^[1 + H2]d 



r2 + r2 



1 



- d^G\ 



de^ +r^ sin 



£> .9G 
K + cot9- 



(142) 
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The Einstein field equations (|45() give the following constraints on the undetermined angular functions ki{9), ^2(6') 
and ^3(6'): 



ki{9) — ciCos9 + C2 
hW + ksie) = C3 cos 61 + C4 



1 + - cos 6 In 



■ cos 9 In 





— cos 9 \ 


ir 


+ cos 9 J 




— cos 9 ^ 




+ cos 6* / 



cos 9\n{l — cos 9) 



(143) 



where Ci, i — 1, . . . ,4 are arbitrary integration constants. Making the choice C2 = = C4 and fc3(0) — brings the 
metric (|142ll in the form 



since in this case 



e^^ll - HQ\dt^ + e-''=[l + iJzldr^ + [l + iJz] d6|2 + sin^^ [l + Hq] 



d^G - - dG 

X + COte— ^ilo 



(144) 



(145) 



4. Comparison with the Weyl class two-body solution 

The solution (| 144(1 we have found to first order in the perturbation is just the linearization with respect to m of the 
exact solution given in Appendix A representing the superposition of a Chazy-Curzon particle and a Schwarzschild 
black hole or of two coUinear Schwarzschild black holes (to linear order these two solutions agree) 



dS^ = -e"^ [1 ~ K^]dt^ + e-"^ [1 + }i^]dr^ + [1 + h'^]d9'^ + ; 



(146) 



where 



hX = 
K = 

K = 



2m 



Zw 



AMm 



1 - 



M — b cos 9 



(147) 



and 



fs{r) = 1 - 



2M 



= [{r-MY + -2{r-M)b COS 9- sin' 9] 



2.11/2 



In fact the metric 1(144(1 has the same form of the Weyl metric ((146(1 : direct comparison between Eqs. p34(l and 1(136(1 
and Eq. ((147(1 shows that the functions Ho and H2 coincide exactly with the corresponding ones Jiq and after 
shifting the location of the particle to b — > b — M. and recalling the definition of the "active gravitational mass" of 
the particle by ,3^J resulting in the overall factor /s(&)^^^ (so that to — > to/s(6)~^/^). 

It is worth noting again that this solution is characterized by the presence of a conical singularity on the polar axis 
between the bodies Hu Hg . 



5. Relation with the Regge-Wheeler approach 

Using Eq. 1(52(1 expressing the metric perturbation functions in the Regge-Wheeler gauge as combinations of metric 
perturbations expressed in an arbitrary gauge together with the conditions 1(100(1 identifying our gauge choice yields 
the explicit relation between the two gauges 



H, 



(RW) 



(BGR) 



ijf = iff + rir - 2A^)G(bg^)" + (2r - 3M)G(bg^)' , 

^(RW) ^ ^(BGR) _^ (y, _ 2M)G'^^^^'>' . 



(148) 
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Since the parameter / does not appear explicitly, the previous relations remain valid for the corresponding summed 
quantities as well, which will be denoted by a bar. The perturber metric written in the Regge- Wheeler gauge in then 
given by 



sm( 



(149) 



where 

where Hq and K are given by Eqs. p34|) and p41|) respectively and Eq. p39|l has been used. 



(150) 



V. PERTURBATION ANALYSIS: THE REISSNER-NORDSTROM CASE 



A. The general perturbation equations 

We now turn our attention to the charged Reissner-Nordstrom black hole case. The geometrical perturbations ft.^^ 
for the electric multipoles are given by 



IVII = 



e'^HoYio HiYio ho 
sym e^''H2Yio hi 



dVi 



10 



89 
80 



sym 



sym 



sym r2 KYio + G 



80^ 



sym 



sym 







sin^ I KYio + G cot 



,8Y, 



10 



39 



(151) 



where the symbol "sym" indicates that the missing components of h^^i, are to be found from the symmetry h^j^^ = 
and e'^ = f{r) is Zerilli's notation. 

In order to simplify the description of the perturbation, we use the Regge- Wheeler j^] gauge to set 



thus we have 



ho = hi = G = ; 

e-'HoYiQ HiYio 
HiYio e--'H2Yio 
r'KYio 



r^sin^OKYio 



(152) 



The electromagnetic field harmonics /^^ for the electric multipoles are given by 





foiYiQ fo2 



961 

antisym ^ 

antisym antisym 

antisym antisym antisym 



(153) 



where /^^ denotes the 0- independent part of /^i^, and the symbol "antisym" indicates components obtainable by 
antisymmetry. 
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The expansion of the source terms (|17|) gives the relations 



(154) 



with 



1 qVW+l 



Sir-b) , 



where we have used the foUowing expansion for S{cos9 — 1): 

S{cose - 1) = Vtt^ V2rTTYio 



(155) 



(156) 



The independent first order perturbations of the quantities appearing in the Einstein-MaxweU field equations 
are given by 



Goo 
Gil 

G12 
Goi 

Gq2 



1 



2 

26" 



2K" -IH2' + -]k' -2(^ + -]{Ho + H2) 



r \ r J 

[{X + l)H2-Ho + XK]^Yia , 
1 r2 / 2\ 2p^'^ 



2 



ly' 1 



+ + + ]{K- H2) \Ym + -iHo-H: 



2 r 
1 



-e'^ { K" +{ ly' + -] K' - Ho" -{- + -]H2' -{— + -] Ho' + 2(A + l)-^{Ho - H, 



Hi\yio , 









^H2- 













Ho 



+ XK] 


'Yio , 






-(¥ 




d^Yio 








dYio 




' 80 ' 





A / , 1 

r 

dYi 



V r 
' I . / 



y{ffi +.H,}—, 



(157) 
(158) 

(159) 
(160) 
(161) 
(162) 



Too — — -r— 
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rripart 
-'00 



F' 



1 Q'e''H2 , .QeVoi 



Stt 



o JOl > J-IO 



Toi = -i^{2^iTiK., 



77— ^oo^;o , 

2~ /(/ + l)e-/o2 

— Joi + r/01 2 



^5{/l2' + .7l2} 



91io 
90 



1 jQ2e--i7o Oe-^foi 



1 J ,Qe-Vo2l 5II0 



Stt 



^K' } Yo , 



* 77i3l>' 



7o2 - -J— S 2— /12 

I {I + 1)6" 



86 ' 



no 



96* ' 



2 /l2j Yio , 



-T— ^ |/oi - /02'| 



>lo ,(163) 

(164) 

(165) 
(166) 
(167) 
(168) 



where X ^ ^ {I - 1) {I + 2) and a prime denotes differentiation with respect to r. 

The angular factors containing derivatives vanish for I = 0; moreover, the two angular factors in the expression 
(|159|l for G22 are not independent when Z 1 (in fact, 8^Yio/89^ = —Yio)- Therefore, the cases Z = 0, 1 will be 
treated separately. 
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For all higher values of I, the Einstein-Maxwell field equations (O imply that the corresponding curly bracketed 
factors on the left and right hand sides are equal, so that the system of radial equations we have to solve is the 
following: 



2K" - - W' 



^ tw' + -\K' - 



= K" 



6 

+ - 

r 



K' 



{W -K)-2 



W 



-W' + K' ~ v'W + 4 



-]K' -W" -2[v' + -]W' 
Qe-Vo2 



JQI + -JOl 



Q , 

^2 



/oi ^ /02' , 



+ 4- 



01 



2v' 

— ]{K~W)-2 



(170) 



foi ,(171) 
(172) 

(173) 
(174) 



since 



Hn — Ho 



W , 



Hi = 



/l2 = 



(175) 



We are dealing with a system of 6 coupled ordinary differential equations for 4 unknown functions: K, W, /oi and 
/o2- Compatibility of the system requires that these equations not be independent. We see that the Eq. 1)1711) can be 
obtained from Eqs. (|170|1 and 1] 172(1 : in fact it is enough to consider Eqs. I|170l) and (|172|l . solving them for K' and 
W, together with the corresponding equations obtained by differentiation with respect to r, solving those for K" and 
W" and finally substituting all these quantities into Eq. ((1711) . which is then identically satisfied. 

Next another equation must be eliminated. It is useful to introduce the the new combinations 



X = K-W , Y ^K + W , 
so that the system we have to solve for the unknown functions X, Y, /oi and /02 is the following 



(176) 



= 

= 




X" + Y" -^{X' -Y')+{^ + ^] (X' + Y')-2 
— H /loo , 

2Ae- 



1 



(X-Y) 



2\p^ d^e^ 
'-^Y-^{X-Y) 



-Ix'-Ux' + Y') 
r 2 



-X - 



01 



X'-i-{X-Y)^S^ 
2 

~ 2~ l{l + l)e~'^fo2 

JQI + -JOl ^ 

/oi — /02' • 



§,{X' + Y') 



Attv 



(177) 

(178) 

(179) 

(180) 
(181) 



First we solve algebraically Eqs. ((178(1 and ((179(1 for /oi and /02 respectively. Then we substitute the quantities so 
obtained into Eqs. I(177() . I(180() and p81() . involving only the gravitational perturbation functions X and Y 



foi 

fo2 







2r2 - 3Mr + _ q2 2Xr^ + 

— A H ■ — — Y — -— — A 



4rg 
- 2Mr + 



4rQ 



Q_ 

4^2 



Y 



4Q 



Y" + -Y' 

r 



2(A + 1) 



4rQ 



:Y 



loo 



= Y 



2 - 2Mr + g2 
(r^ - 2Mr + Q^)X" + 4(r - M)X' - 2XX , 
„ , 2 2(A + 1) , IGirrQ 



2{r + Q)ir~Q) 

(r2 - 2>lr + 02)2 r{r^ - 2Mr + Q'^) 



X' 



:Y 



2{r + Q)ir-Q) 



X' 



2(2A + 1) 

^ r2 - 2Mr + 

2(2A + 1) 



r-2-27Wr + 02 {Mr - Q^) r{r^ - 2Mr + Q^) - 2Mr + 



(182) 

(183) 

X ,(184) 
(185) 
(186) 
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B. The compatibility and Bonnor's condition 

Direct comparison between Eqs. H184fl and H186|) shows that they are compatible only if the following stability 
condition holds 



The same condition can be obtained by imposing conservation of the stress-energy tensor: 



= Aoo - 327r^z; . (188) 



So Eq. (|169|l becomes 



= 6^'' 



2K" ^lw'+(iy' + -]K'-4(X + -\W 



(^W + K)-2'±-, 4^1-i2i-327r5-.. (189) 



This equation is identically satisfied, seen by substituting into it the quantities K' , W' , and K" and W" obtained 
from Eqs. (|170|l . (|172|l and their derivatives, and /oi' and /02' from Eqs. (|173|l and (|174|l respectively. 

The stability condition (|187|l coincides exactly with the equilibrium condition H44|l for such a system, which has been 
discussed by Bonnor 19] in the case of a test field approximation. Therefore the static configuration of a particle at 
rest near the black hole remains an equilibrium configuration as a result of the perturbation only for certain positions 
of the particle which are completely determined by the charge to mass ratio of the black hole and of the particle itself: 
in particular, if we require that Q / M. < 1, the particle must be overcritically charged, i.e. q/m > 1. However, the 
choice of equilibrium configurations is a special one, since it forces us to consider particle and black hole charges of 
the same sign only. Hence, if we want to consider the more general case of charges of opposite sign, the nonvanishing 
contribution to the stress-energy tensor due to some external force should be taken into account in solving the full 
Einstein-Maxwell system (|SJ. In the extreme black hole case {Q — Ai), instead the equilibrium condition H187|) 
becomes simply 

m^q , (190) 

implying that all the configurations arc equilibrium configurations as a result of the perturbation, but the particle 
also is forced to have the critically charged ratio q/m — 1. 

1. Discussion of the equilibrium condition 

Let us examine the equilibrium condition (|187l) in more detail; it can be rewritten in the form 

q Q b/M-iQ/M)^ 



mM [(6/X)2-26/7W + (Q/7W)2]i/2 



T{b/M;Q/M), (191) 



with T{hl M.\QI M) > 1 for all values of jQ/A^j < 1 and b > r+. Thus the system turns out to be coupled in such 
a way that, for a fixed value of Q/M, the charge to mass ratio q/m of the particle must decrease as the separation 
parameter b/A4 increases, while it must be larger and larger as the particle approaches the black hole horizon. On 
the other hand, for a fixed distance parameter b/M, the value of q/m must increase for decreasing values of Q/M, in 
order to oppose the gravitational attraction and so maintain the equilibrium, while it must decrease until it reaches 
the value 1 as the black hole approaches the extreme condition, since J-{b/M;Q/M ^ 1) ^ 1 irrespective of the 
value of b/M. Obviously this is the limiting situation of two critically charged bodies which are in equilibrium at any 
separation distance. 

Summarizing (see also Fig. ^: 

a) Q/M fixed: 

for 6 — r_(_ , q/m — > 00 ; 

for 6^00, q/m^M/Q; (192) 
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b/M Q/M 



(a) (b) 

FIG. 1: Fig. (a) shows the behaviour of the charge to mass ratio q/m of the particle as a function of the separation parameter 
h/M, according to the equilibrium condition 119111 . for a fixed value of Q/M = 0.6. For 6 — > oo we have that q/m — > M/Q ~ 
1.67. Instead in Fig. (b) the parameter q/m is plotted as a function of the charge to mass ratio Q/M of the black hole, at a 
fixed distance b — 8M. For Q/M — > 1, we have that q/m — + 1 as well. 

b) b/M fixed: 

for Q/M~^0, g/m— >oo; 

for Q/M ^ 1 , q/m 1 . (193) 



C. Solutions of the perturbation equations: the extreme case {Q = M) 

Consider first the cases I — 0, 1. 



1. The I — case 



The relevant equations come from the first order perturbation quantities (|157() - (|168|l appearing in the field equations, 
which do not contain angular derivatives and which reduce to the foUwing for I — 

r r 7-2 2{r - My 

-2M~fm , (194) 

= ^ ^r^Ho' -{r- M)K' + H2-K- —Ho + 2Mfm , (195) 

= (r - Mf[K" - Ho"] + {r- M)[2K' - H^'] - ^ Ho' - ^i?2 + 4M/oi , (196) 

= /oi' + ^/oi - ^K' + Anv . (197) 
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Let us look for a solution of the form Hq = H2 = W and K = W . Then the preceding equations reduce to the 
following 



= ir^MfW"+'''-'^:^^\' + ^W- 



2{r-M)' 



01 



-W - 2M/01 , 



= /oi +-/o 



W + 47rw 



By solving Eq. p99|) for /oi and substituting it into the other equations we obtain 



= W" + -W' 



^4 



r 2{r-MY 

Q ^ W" + -W + . 

r r — M 



^00 , 



These are actually the same equation, as follows from Eqs. (|155|l and H190|) . and the solution is simply 



W = 4\/^- 



r - M , b- M ^, , ' 

- r) H {}{r - b) 

r r 



b-M 

after a suitable choice of the integration constants; then from Eq. p99|l we obtain 



foi = 



r — A4 

2M- — —d{b - r) - (r - 2M)d{r - b) 



(198) 

(199) 
(200) 

(201) 
(202) 

(203) 
(204) 



It is worth noting that the function /02 turns out to be undetermined since it does not appear in the system (|194f) - 
(|197|l . However, the tensor harmonic expansion H153|) of the electromagnetic field shows that the angular factor 
corresponding to the component /02 vanishes identically for I = 0. 



2. The I = 1 case 



The derivation of the relevant equations from Eqs. I|157|l - (|158|l requires particular acare in this case. In fact, the 
two angular factors in the expression (|159|l for G22 are not independent when 1 = 1, since d'^Yio/dO'^ = — Yio, so the 
two separate terms in the corresponding field equations collapse to one. The relevant equations are thus given by 







2(r-X)2 
(r-X)2 



^00 - 2Al/oi , 



Ho' - {r - M)K' + H2 



1 



= {r-Mf[K"-Ho"] + {r-M)[2K'-H2'] 



Ho + 2Mfoi 



7-2 + Mr - 2M'^ , 
Ho + Ho 



1 + 2- 



H2 + 'iMfoi 



= {r-Mf[K' -Ho']-ir-M)H2 



- 3Mr + 2M'^ 



Ho + ^Mf, 



02 



— /oi' + -/oi — , . .s 
r [r - M) 

= /oi — /02' • 



2/02 -^K' + At:v 



(205) 
(206) 

(207) 

(208) 
(209) 
(210) 
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We look for a solution of the form Ho = H2 = W and K = W. Thus the preceding equations reduce to 



, 2r'-3Mr + M^^, 



{r - MyW" + 



2- 



W ^ 



2{r-My 



:Aoo - 2Mfoi 



r - M 



W + 2f, 



02 



2 ~ 2 
= /oi' + -/oi- (^—^ 



— /oi ^ /02' ■ 

By solving Eq. H212|) for /oi and Eq. (|213|l for /02, and substituting the results into the other equations we obtain 



(211) 

(212) 
(213) 
(214) 
(215) 



= W" + -W'- 



= w" + -w'- 



{r-Mf 2(r-MY 
W ^ —V . 



00 



[r-M) 



M 



These are actually the same equation, as follows from Eqs. (|155ll and (|190|l . and the solution is simply 



W = 



4^/7? TO 



r- M\ , b- M 

i9(&-r)H — i?(r-6) 



r-M 



after a suitable choice of the integration constants; then from Eqs. H212|l and H213() we obtain 



Joi 



2Vnq_ 
V3 



r-M 



b-M 



(r + 2M)'d{b - r) - 2(5 - 7W)i?(r - b) 



~ _ 2y/^ r-M 



r — M \ b — M 

^ z9(6-r) + — ^i9(r-6) 



b-M 



r-M 



(216) 
(217) 

(218) 

(219) 
(220) 



Thel> 2 case 



Finally, consider the general case I > 2. The system of equations to be solved for the unknown functions X, Y, /oi 
and /02 coming from Eqs. 1178(l - (|181|l is 



~ _ 2r^ - 3Mr + M^ , r-M , 2\r^ + M^ M 
^ AM? AMr^ 



/( 



02 



AM " ' 4r 
= {r- MfX" + A{r- M)X' -2XX 



2(A+1) 



Y 



Idtirr 



= F" + -Y' . , 

r (r — MY {r-M) r{r - M) 



2{r + M)^, , 2(2A + 1) 
{r-MY 



X 



where the relation (|190|l must be taken into account once the solution is found. 
Eq. H223(l involves the function X only and can be solved exactly 



X ^ ci{r - My-' + C2{r - M) 



-{1+2) 



which can be re-expressed as 



X^{r- M)-'[ciTi{r) + C2T _i_i{r)] 



by introducing the following notation 



Ti{r) = {r-M)' 



X 



(221) 

(222) 
(223) 
(224) 



(225) 
(226) 
(227) 
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At this point, it is enough to solve Eq. (|224() to obtain the complete solution of the system. The general solution of 
the homogeneous equation is simply 

r — A4 

11>om = [c3Tz(r)+C4T_,_i(r)] , (228) 

r 

so that the solution of Eq. (|224|1 becomes 

Y = Fhom + B^—^ [Tz(6)T_z_i(r) - T,(r)T„,_i(fe)] i9(r ~ b) + c^niir) + c^n^i-iir) , (229) 
r 

where the quantities B and Qi (r) stand for 

80F 



B = 



V2r+T 

^lir) = ^[X(/2-; + l)T,_i(r) + /T;(r)] . (230) 

The functions f2;(r) and Q-i-i(r) are not regular at the horizon nor at infinity, so that we must set the constants 
Ci = — 02- This means that there is only one solution which satisfies the suitable regularity condition: the solution 
with X = 0, and thus K = W = Y/2, from the relations l(T7^ . 

The remaining undetermined constants C3 and C4 appearing in the solution 1229|l can again be found easily by 
imposing regularity conditions on the horizon and at infinity, giving the following final form for the solution 

Y = b'—^ [Ti{r)T_i_,{b)d{b ~r) + Ti{b)T _i_^{r)d{r - b)] . (231) 
r 

The electromagnetic perturbation functions /oi and /02 are then easily determined from the relations (|221|l and 
(|222|l respectively, which for X = become 

r — A4 

fo2 = ■ (233) 



4- The analytic solution summed over all values of I 

In order to reconstruct the solution for gravitational and electromagnetic perturbation functions for all values of Z, it 
is important to note that the both solutions (|203|) and H218(l for W correspond to the expression H231|l for the function 
Y/2 evaluated for / = and / = 1 respectively. The same holds for the electromagnetic perturbation functions /oi 
and /02: if we set / = 0, 1 in H232|) and (|233|l . we obtain exactly H2U4|) . (|219|l and (|220|) . Therefore, we can take the 
expressions (|231f) . (|232|) and H233|) as the solution for all values of /. 

The gravitational field harmonics /i^j^ reduce to 



e-'^Yio 











e- 






























r2 sin^ e^Yio _ 



(234) 



since we have K — W — Y/2. Now, if we denote the sum of the radial function Y over the spherical harmonics by y, 
it is easy to show that this quantity can be summed exactly using the properties of the Legendre polynomials Pn{x), 
whose a generating function is just 

00 

5(^,2;) = [1 - 2xi + t^]-^/2 = ^F„(x)t" , \t\<l. (235) 

11=0 

In fact, setting t — {b — A4) / {r — A4) and x — cos 9 in the generating function H235|l leads to 

1 = V P„ (cos 9) Z^',-^^" , (236) 

n— 



27 



for r > h, where 



V=[{r- Mf + {b- Mf - 2(r - M){b - M) cos 6] 



1/2 



alternatively, choosing t — {r ~ Ai)/{h — M) implies 



cost 



{r-MY 

n— ^ ^ 

for r < b. Therefore, the fohowing representation formula holds 

-=y 



?i=0 



leading to the result 



I 



l=Q 



{h-My 



P;(C0S( 



r-M 1 



r V 



(237) 



(238) 



(239) 



(240) 



At this point, after recalling the definitions (|176l) for K and W (from which K = W = Y/2, for X ^ 0) and the 
relation we have that the solution summed over the harmonics for the perturbed gravitational field turns out 

to be completely determined by the function 



y r-M 1 

ri = — = 2m — 



2 r V ' 

so that the new line element ds^ from the first of relations and Eq. (|234(l is then 

ds^ = -[l-H]f{r)dt^ + [l + n][f{r)-'^dr^ +r^{de^ +s\t\^ 6d(l)'^y\ 



(241) 



(242) 



It can be shown that this perturbed metric is spatially conformally fiat; moreover, the solution remains valid as long 
as the condition \Ti\ <C 1 is satisfied. 

Next by using the relations H232(l and (|233|l . the knowledge of the gravitational perturbation function Ti. (or y) in 
closed form permits us to give an exact reconstruction of the perturbed electromagnetic field /^^ as well, namely 



f f V r-M dy M 

foi = 2^ foiYio = — + —y 



4r dr Ar^ ' 
~ dYio _r-Mdy 

02- 



de 



4r 89 ' 



(243) 
(244) 



so that the electric field components and Eg are given by 



Eff 



loi 



q{r - M) 



-2M 



r{r - M)[ir - M) - (b - M) cos9] 



X,2 



(245) 



the total electromagnetic field to first order in the perturbations is then (from the second of relations @) just 

'M 



Er 



dt/\dr- Egdt A dO 



Let us verify that Gauss's theorem is satisfied by integrating the dual of the electromagnetic form 
5 containing both the charges Q — A4 and q 



'F AdS = An{Q + q) 



(246) 
over a surface 

(247) 
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We only need to calculate the component *Fe<f>, which results to be (to first order in the perturbation) 



sin ( 



(248) 



Let us calculate separately the two different contributions to the integral (|247|l due to the background and particle 
electric fields. The former one is given by 



$RN = 27r / r'^ sm9{l +n) — de ^ AnM + ArrMq / 



sin0 



-de 



(249) 



recalling the definition (|241|l for the gravitational perturbation function H. Next it is easy to show from (|237() that 



V 



-de 



1 



dV 



de : 



X>(7r) -P(0) 



{r-M)ib-M) Jo de (r - M)ib - M) 
^r-M)ib-M) ^^^^^ -r) + ib- M^r - b)] 



Thus the flux l|24H|l becomes 



$RN = 47r7W + -|^^^^^[(r - M)^9ib - r) + {b - M^r - b)] . 



r{b-M) 

The charged particle contribution to the flux is given by 



$pa^rt = 2tt I r'^ sine Erde ^ -Anr'^-^ 



(250) 



(251) 



r - M 



4r 



T-Lsinede 



d_ 

dr 



smt 



■de 



so that using the expression (|250|l it reduces to 



'inq — 



r — M 



Therefore the total flux $ is 



* = *RN + ^'part = ^T^M + 47rgz9(r - b) 



(252) 



(253) 



(254) 



so Gauss's theorem (|247f) is satisfied. 

As a matter of fact there exists an alternative way to treat the problem, at least for the case of particle and black 
hole charges of the same sign, which permits one to deal ab initio with closed form expressions both for the perturbed 
gravitational and electromagnetic field, as we will see in the next section. That leads to enormous simplification in 
the treatment of the perturbations of extreme black holes, and it is a powerful tool for the analysis of higher order 
perturbations, since one can stop in the perturbation theory at any order one desires. 



D. Perturbations of the extreme black hole: an alternative approach 



Although the spacetime of a black hole and a charged particle is not an electro-vacuum one, we start from an exact 
solution of the Einstein-Maxwell field equations without sources, and by linearizing it, we will obtain an exact solution 
of the coupled perturbation equations. 

The exact solution from which we start is the Majumdar-Papapetrou solution describing a system of many extreme 
Reissner-Nordstrom black holes. We note that by using the Israel- Wilson-Perjes approach one could also include 
rotation in the present formulation. In particular, we will be interested in the two body problem, which is discussed 
in detail in the monograph of Chandrasekhar j33] . 

The physical scenario we have in mind is the one in which one black hole has a mass (charge) absolutely smaller than 
the one of the companion object. In this limit and assuming a separation distance of the two objects not comparable 
with their gravitational radii, we can interpret the system as an extreme charged black hole interacting nonlinearly 
with a very small charged object. The charge-to-mass ratio q/m = 1 of such a particle is not a problem for the 
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effect we are going to discuss. In fact in the equations of Zerilli the charge and mass of the particle appear in the 
source terms of the coupled ODE system only as a common multiplicative factor. The particular integral built with 
the solutions of the homogeneous system will be proportional to such quantities, and this makes the mathematics of 
the problem (but not the physics, of course) naively insensitive to a particular choice of such parameters, provided 
always that the perturbative regime is valid. Clearly this assumption is very delicate, because the black hole with a 
particle is not a vacuum spacetime and we mimic it with an electro-vacuum one. In some sense the idea is to use an 
analogy such that the particle and a very small extreme black hole can be thought as the same kind of object and 
have the same fields at far distances from themselves. The Majumdar-Papapetrou solution is included in the class of 
electro- vacuum and conformally static solutions of the Einstein-Maxwell system. 

The three-dimensional conformally flat geometry in this case is allowed only for extreme black holes, and this fact 
demonstrates the special nature of Reissner- Nordstrom black holes with Q ~ Ai. In this situation the complicate 
system of non linear PDEs collapses into a simple complex Laplace equation, which acts as master equation for 
building up infinite solutions; however, apart from the Majumdar-Papapetrou spacetime, the other solutions give rise 
to naked singularities. It is important to note that, because of the linearity of the Laplace equation, a superposition 
of many relativistic bodies is possible. 

We are now ready to discuss quantitatively the Majumdar-Papapetrou spacetime. In fact, by introducing the 
quantity 



M 



it is easy to show that a solution of the electro- vacuum Einstein-Maxwell system is given by 



U'^[dR^ + R^{d0^ + sin^ ed(t>^)] 



(255) 



(256) 



where U satisfies the 3-dimensional Laplace equation V^J7 = 0, and the electrostatic potential is given by = U'^. 
The quantity b is the separation distance between the two bodies, and the quantities M, m can be identified with 
the masses, if the sources are far enough apart. The charges of the system are equal to these masses and have the 
same sign. This is a limit of the analogy we are using in this section, because we can study an extreme black hole and 
a charged particle provided that both have the same sign of the charge. Bodies with different signed charge would 
require the study of the full Zerilli system, as we have done in the previous section. 

The manifold under consideration has a coordinate interpretation which is absolutely nontrivial, in particular in 
relation with the location of the horizons and the singularities. However, because we will be interested in the study of 
the perturbations of such an object, we can use a coordinate transformation R — r~A4, which in the case m = gives 
immediately the usual Reissner-Nordstrom solution, and consequently we can construct the gravitational perturbation 
with respect to m referred to such coordinates by expanding the metric to first order in m: 



,MP _ 



"■rr 



.MP 



2m (r - Mf 

2m r 
E r-M ' 

—r{r~M) , 
r{r — M) sm 



(257) 



where 

S = [(r - Mf + {b- - 2(r ~M){b- M) cos 9] ^'"^ 
and the following representation holds (see Eq. (|239|l ') 



(r) 



20F 

V2rrT 



moreover, recalling the notation l|227|l introduced in the previous section, 



one finds 



Ci{r) 



V2r+T 



[T,(r)T_,_i(6)z9(6 - v) + Ti{h)T _i.^{r)'d{r - h)\ 



(258) 



(259) 



(260) 



30 



We can now suitably scale the gravitational perturbation functions (|257|l as in the Zerilli approach 



^MP _ j,MP -i/I j,MP 2T, __ uMP 

finding that in this case 



(261) 



with 



- 2m r-M 2m(r-A^) ^ 

/I = -^^ = yCi{r)Yio . (262) 

I 

This quantity is equal to the gravitational perturbation function (|241() obtained in the preceding section following the 
approach of Zerilli, since g = m in the present treatment, and !] = !?. 

The nonvanishing components of the effective electromagnetic tensor (that is, the total electromagnetic tensor, from 
which the contribution Ffi^^ = —{M./r'^) dt A dr of the "background" field of the black hole of mass M. is subtracted) 
are easily evaluated from the metric (|256|l and the vector potential = U^^dt 



p 



MP 



2r 



m{r - M) E[r(r - 4A^)E - 2mAl(r - M)] + r^[ir - My ~ {b ~ M) 



MP mjr^M) jb-M) 
^fo — ft; 7 rTTT^smt^. (263) 

By expanding the previous expressions to first order in m we also obtain the electromagnetic perturbation 

MP _ m{r-M) \ , r{r-M)\r-M-{b-M)cos.d\ 

Jtr - 1 ^-^^l ^2 



/MP = (b-M) Sine. (264) 



m(r — A^)'' 

As in the case of the metric perturbations, the previous quantities are the same as those (see Eq. ()245|l ^ obtained in 
the preceding section, with the identifications 

Although the metric perturbation /ijfjf continues to be singular for r — > as in the background case, inspection of 
the curvature invariant K. = RapR"^ shows that its perturbative expansion with respect to m there gives 
and /C*-^-* = 0. Moreover, the electromagnetic perturbation fjfj' in the same limit is zero. This fact, in agreement 
with the "No-Hair Theorem", implies that the extreme black hole does not feel the external charged particle and 
preserves the shape of its horizon. This is clear from the stability theorems for the electro-vacuum perturbations of 
the Reissner-Nordstrom background, which require regularity on the horizon for all the modes and consequently for 
the fields built from them. 

We have thus demonstrated the complete equivalence between this new approach with the standard one of Zerilli; 
in addition, the result is more general, since the Majumdar-Papapetrou solution is an exact solution of the Einstein- 
Maxwell equations. 

It is possible to reach the same conclusion directly, by verifying that the quantities H257I) and H264|l satisfy the Zerilli 
equations for the Reissner-Nordstrom perturbations, as shown in the following. The perturbed effective electrostatic 
potential is easily determined, by expanding Af^^ = to first order in m, and by subtracting the contribution 

(1) — -M/r due to the black hole of mass Ai: 

ar = 1 - ^^^^^ = 1 - ^^^^^ E ^^ir)Y,o , (265) 

I 

where the latter expression is obtained by expanding the former one in spherical harmonics. Consequently the 
electromagnetic perturbation function is given by 

/MP ^ d_^MP (266) 

At this point, denoting by ftr and h the 6'-independent parts of the expansions 1)2661) and (|262|) respectively (so that 
the meaning of the notation fte is also clear), it is easy to show that these quantities satisfy the Zerilli equations 
(|170|l . H172|l . H173|) and (|174|l for the extreme case together with the stability condition H190|) . with the identifications 



q = m , K = W = h , f^^ = f^.^ ^ /o2 = /* 



te 
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E. Solutions of the perturbation equations: the general non-extreme case (Q < A4) 

Let us consider first the cases / = 0, 1. 

1. The I = case 

The relevant equations come from quantities (jl57|l - Hlfi8|l which do not contain anguiar derivatives, and which in 
this case become 



2{r^-2Mr + Q^) 



Aoo - 2Q/oi , (267) 



= \: Ho' -{r- M)K' + H2-K- ^Ho + 2Q/oi , (268) 

= (r2 - 2Mr + Q^)[K" - Ho"] + (r - M)[2K' - H^'] - !:1±^±_Z^Ho' 

r 

002 

-^i?2 + 4Q/oi, (269) 

= foi + ^/oi - ^K' + inv . (270) 

Let us iook for a solution of the form Hq = H2 = W and K — W . Then the preceding equations reduce to the 
following ones 

= (r^ - 2Mr + Q')W" + ~ ^"^'^ + ^ V ' + ^VF + . , ^ Aoo 

r H 2{r^ — 2Mr + Q'^) 

-2g/oi , (271) 

= \, + - 2Q/01 , (272) 

= hi + Ifoi - + Attv . (273) 

By solving Eq. I|272|l for /oi, and substituting it into the other equations we obtain 



= W" + -W + -:P^v . (275) 

r Mr - 



These are actually the same equation, as follows from relations (|155|l and (|187|l . and the solution is simply 

W = iy/^^f(b)-^/^- Ur - M)^(b -r) + (b- M)d(r - 6)1 , (276) 
r 

after a suitable choice of the integration constants; moreover, from Eq. I|272|l we obtain 

/oi = 2V^ ^^'^^^ -^{ [M{Mr - 2Q2) + ^g2] ^(5 _ ^) _ _ 2Q^){h - M)^{r - 6)} . (277) 
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2. The I = 1 case 

The relevant equations coming from quantities H157() - (|168() are given by 

.2^^w/ 3r2-57Wr + 2g2 - 2Mr + 



= [r^ -2Mr + Q^)K'' 



Q2 
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/02 - + 47ru 
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(279) 

(280) 

(281) 
(282) 
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Let us look for a solution of the form Ho = H2 ^ W and K ~ W. Thus the preceding equations reduce to the 
following ones 



= {r^-2Mr + Q^)W"+ ^'^^ ^"^^ + ^ V ' 



2 - ^ 



W + 



2(r2 -2A^r + g2 



-An 
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(284) 
(285) 

(286) 

(287) 
(288) 



— /oi — /02' ■ 

By solving Eq. (|285|l for foi and Eq. H28t)|) for /02, and then substituting the results into the other equations we obtain 



= W" + -W' 



= W" + -W'- 



rW 



-2Mr + Q^ 2{r^-2Mr + Q^y 
tW + — —V 
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(289) 
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These are actually the same equation, as follows from relations (|155|l and (|187|l . and the solution is simply 
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(291) 
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after a suitable choice of the integration constants; then from Eqs. 1285|l and H286() we obtain 

%/3^ q b-r^ r Mr{r^ - 3Q^) - 2rlrl 
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(292) 



3. Thel>2 case 

Let us look for the solution of the system consisting of Eqs. I|182(l . 1)183(1 . 1(185(1 and 1(186(1 for the general non-extreme 
case, together with the relation p87l) . 

Eq. ((185(1 involves only the function X, and can be solved exactly. Putting X — f{r)^^^^w{r)/r and making the 
transformation z = {r — A4) /T Eq. ((18511 becomes 



= (1 - z^)w" - 2zw' + 



l{l + l) 



1 



1 



(293) 



where primes now denote differentiation with respect to the new variable z. The general solutions are the associated 
Legendre functions of the first and second kind Pi'{z) and Q}{z), implying that 



X 



fir)-' 



[cifi{r) + C2gi{r)] = CiXi + C2X2 



(294) 



where the functions fi{r) and gi{r) have been defined in 1(38(1. The arbitrary constants ci and C2 must be set both 
equal to zero in order that the function X to satisfy regularity conditions on the horizon and at infinity, in agreement 
with the analysis of the solution we have done in the extreme case. At this point, it is enough to solve Eq. ((185(1 with 
the condition X = and thus K = W = Y/2 to obtain the complete solution of the system. 
Let us consider Eq. I(186() in which we set X ~ 0, namely 



= (r^ - 2Mr + Q^)Y" + -(r^ - 2Mr + Q^)Y' - l{l + \)Y 



lQ'Kb^f{b)Q 
Mb-Q"^ ' 



Putting Y — f{rY/'^w{r) and making the transformation z = (r — A^)/r the previous equation becomes 
= (1 - z'^)w" - 2zw' + 



S{z - fi) 



(295) 



(296) 



where primes now denote differentiation with respect to the new variable z and (3 ^ [b — M)/T. This equation is 
formally identical to Eq. I(31() apart from a different coefficient in front of the delta function. The general solutions 
of the corresponding homogeneous equation are thus the associated Legendre functions of the first and second kind 
PI{z) and Q\{z). Taking the functions l(38() as the two linearly independent solutions of the homogeneous equation, 
the solution of Eq. ((296(1 turns out to be given by 



Y = -oVTi'TTi— rv9 
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Note that for Z = 1 this solution reduces to Eq. H291() . since W = Y/2. 

The electromagnetic perturbation functions /oi and /02 are then easily evaluated from the relations (|182|) and H183|) 
respectively, which for X = become 
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4- The analytic solution summed over all values of I 

As in the extreme case we are able to completely reconstruct the solution for all gravitational and electromagnetic 
perturbation functions in closed form also in this case. 

Consider first the gravitational perturbation function Y, whose solution is given by Eq. (|297|l for I > 1 and Eq. 
(|276|l . with W = Y/2, for I = 0. The sum over all multipoles turns out to be 
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where the quantity T) = Z'rn is equal to 
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and the representation formula 134|l has been used. Since K = W = Y/2 and taking into account the relation (|187|l . 
we have that the solution summed over the harmonics for the perturbed gravitational field turns out to be completely 
determined by the function 



so that the new line element ds^ from the first of relations and Eq. (|234|l is then 

dS^ ^ -[1 - H]f{r)dt^ + [1 + n][f{ry'^dr^ + r^{de^ + sin^ ^d^^)] 



(302) 



(303) 



The asymptotic mass measured at large distances by the Schwarzschild-like behaviour of the metric of the whole 
system consisting of black hole and particle is given by 



Meff ^ M+m + E; 



where the interaction energy turns out to be 



M 



int 



f{h) 



-1/2 



(304) 



(305) 



It can be shown that this perturbed metric is spatially conformally flat; moreover, the solution remains valid as long 
as the condition |7Y| ^ 1 is satisfied. Let us choose the mass m itself as a smallness indicator. FigureElthus shows 
the regions where our perturbative treatment fails, assuming that \TL\ ~ m/M. 
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Introduce isotropic coordinates {t, p, 0, (f>} such that 



r(p) = M+ p + 



The perturbed metric H3U3|) then becomes 

dS^ - -[1 - ■Hip)]A{p)dt'^ + [1 + + p^{dO^ + sin^ 6id02)] 



where 
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By inspection of (|3U7fl we see immediately that the spatial 3-mctric 
is conformally flat, with conformal factor 

^=[l + H{p)]B{pf . 

The condition gtt = defines the new coordinate horizons for the system "black hole + particle" : 

= /(r)[l-7Y] , 
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(309) 

(310) 
(311) 
(312) 
(313) 



where the quantity Ti is given by Eq. H302|) . In addition to the usual black hole inner and outer horizons r_ and 
r+ there would exist another one, implicitly defined by the equation = 1, as a function of the variables r and 9. 
Hence, the presence of the black hole seems to induce the formation of a small horizon around the particle, which 
however would lie inside the nonaccessible region (see Fig. |5J)- Therefore, such a conclusion cannot be drawn since 
our (perturbative) treatment is no longer valid there. 

The reconstruction of the perturbed electromagnetic field /^^ by means of relations (|298|l and 1)2991) and by using 
the closed form expression (|302ll for the gravitational perturbation function H (or Eq. H300() for y) is the following 
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so that the electric field components and Eg are given by 
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the total electromagnetic field to first order in the perturbations is then (from the second of relations ^) 
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FIG. 2: The surface \'H\ = m/M is plotted as a curve in the polar coordinate plane with the choice of parameters Q/ M — 0.6 
and b/A4 = 6 (dashed curve). The whole region inside this limiting surface is not accessible since the condition |7i| <C 1 is 
not satisfied there (in fact, in this region we have that |7i| > m/M). The black circle represents the black hole horizon. The 
"induced" horizon around the particle location is also shown (shaded region) for the further choice m/M = 0.1. 

We are left to verify that Gauss's theorem (|247|l is satisfied also in this general non-extreme case. The only non- 
vanishing component * Fg^ of the dual of the electromagnetic form l|317|) is given by (to first order in the perturbation) 



*Fg^ = sin( 



(318) 



Let us calculate separately the two different contributions to the integral (|247|l due to the background and particle 
electric fields. The former one is given by 
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From relations (|187|) . H3U1|) and (|302|) it is easy to show that 
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Thus the flux H319|) becomes 



Relation H315|) leads to an expression for the perturbed electrostatic potential 

Mr-Q' 

where H — y/2 from Eq. H302|l . Hence the contribution to the flux due to the charged particle is given by 
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after using the expression H320|l . Therefore the total flux $ turns out to be 
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V{Tr) - P(0) = 2[(r - - r) + (6 - X)i9(r - b)] , 



(324) 



(325) 



from Eq. H3()l(l . So Gauss's theorem is not satisfied. However, it is enough to add to the perturbed electrostatic 
potential H322|l the following term 
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which vanishes in the extreme case and whose contribution to the flux is just 

r2 



Mb^Q^ 



Hence the total flux <I> is given by 
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and Gauss's theorem (|247|l is satisfied. It is worth noting that the addition of the term H326|l also leads to a change 
in the perturbed metric functions: from Eqs. I|271(l - (|273|) it is easy to show that the solution for W changes simply 
by a constant, implying a modification of the function H by the constant term —2qT^ /[Q{A4b — Q^)]. However this 
term can be eliminated by a suitable gauge transformation of the perturbed metric. 

Let us denote by Vtost the electrostatic potential of the particle alone obtained within the test-field approximation 
by Leaute and Linet |10|. and by 



Vtot = V + V + V 



BH 



(329) 



the total perturbed electrostatic potential obtained by summing (|322|l and (|326() plus the contribution H42|l of the 
black hole itself. Direct comparison between the potentials shows that they are related as follows 
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The second and third terms in the bracketed expression of H330() represent the "gravitationally induced" and "electro- 
magnetically induced" electrostatic potential respectively and the equilibrium condition H44|) has been conveniently 
used. 



F. Comparison with the Weyl class double Reissner-Nordstrom solution 



The solution for two Reissner-Nordstrom black holes belonging to the Weyl class, once linearized with respect to 
the mass m and charge q of one of them, is given by 



where 
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and 



with 
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Direct comparison of this solution with our solution shows that they do not coincide. In addition, the metric l|331l) is 
characterized by the presence of a conical singularity between the bodies, which can be removed only if both of them 
are critically charged (see Appendix B); our solution is instead totally free of singularities. 



VI. CONCLUSIONS 



The problem of the interaction of a massive charged particle at rest with a Reissner-Nordstrom black hole has been 
studied taking into account both electromagnetic and gravitational perturbations on the background fields due to the 
presence of the particle. Following Zerilli's approach to the perturbations of a charged static black hole, we derived the 
corresponding solutions of the linearized Einstein-Maxwell equations for both the electromagnetic and gravitational 
perturbation functions. We were able to exactly reconstruct these functions summing over all multipoles, giving closed 
form expressions for the components of the perturbed metric as well as the electromagnetic field. A detailed analysis 
of the properties of this solution will be presented elsewhere. 
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APPENDIX A: THE WEYL CLASS TWO-BODY SOLUTION 

Axisymmetric static vacuum solutions of the Einstein field equations can be described by the Weyl formalism [36| . 
The line element in coordinates (t, p, z, cj)) is given by 

ds^ = -e'^dt' + e^f^-'^) [dp" + dz'] + p^-^^dcj,' , (Al) 

where the function ip and 7 depend on coordinates p and z only. The vacuum Einstein field equations in Weyl 
coordinates reduce to 

= 1p,pp + -^1p,p + Ip^zz , 

= 7,p - P[^l - ^l] , 

= 7,^ - 2p^,p^^z . (A2) 

The first equation is the 3-dimensional Laplace equation in cylindrical coordinates; so the function plays the role of a 
newtonian potential. The linearity of that equation allows to find explicit solutions representing superpositions of two 
or more axially symmetric bodies. In general, these solutions correspond to configurations which are not gravitationally 
stable because of the occurrence of gravitationally inert singular structures ( "struts" and "membranes" ) that keep 
the bodies apart making the configuration as stable. In the case of collinear distributions of matter displaced along 
the symmetry axis, this fact is revealed by the presence of a conical singularity on the axis, the occurrence of which 
is related to the non-vanishing of the function 7(p, z) on the portion of the axis between the sources or outside them. 
In the former case we can interpret the singular segment of the axis as a strut holding the bodies apart, while in the 
latter one as a pair of cords on which the bodies are suspended. 

For the static axisymmetric vacuum solutions the regularity condition on the axis of symmetry ( "elementary flat- 
ness" ) is given by 

lim 7 = . (A3) 

p^O 
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a. Superposition of a Schwarzschild black hole and a Chazy-Curzon particle 



The solution corresponding to a point particle of mass m at rest on the symmetry axis sX z = b above the horizon 
of a Schwarzschild black hole with mass M at the origin is given by the metric (jAl|) with functions 



^/j = V^s + -iACb , 7 = 7s + 7Cb + 7sCb , 



where 



and 



Ri 



m 

i?Cb ' 



2^2 



7Cb 



1 m p 



2 i?4 

t-'b 



(A4) 



(A5) 



2 



R1+R2- 2M 



R1+R2 + 2M 



7s 
R2 



2 



(Ri + R^f - 



4i?ii?2 



Vp2 + (z + A^)2 



v/p2 + (z-Al)2 , 

while 7sCb can be obtained by solving Einstein's equations (|A2|I 

m {h + M)Ri - [b - M)R2 



7sCb 



-RCb 



C 



(A6) 



(A7) 



The value of arbitrary constant C can be determined by imposing the elementary flatness condition ljA3|l . We can 
make 7sCb zero between the mass and the Schwarzschild source by choosing C = —2mM/{b^ — M^), or outside them 
by taking C to have the opposite sign, but C can not be chosen so that 7sCb vanishes on the whole z axis, giving rise 
to the well known conical singularity, corresponding to a strut in compresson which holds the black hole and particle 
apart. 

The relation with standard Boyer-Lindquist coordinates in the Schwarzschild case is simply 



P 



2Mrsm9 , 



z — (r ~ M) cost 



In Boyer-Lindquist coordinates, under the transformation (|A8|) the metric (jAl|) re- writes as 



^2 = -e^^dt^ + e^^-'-^^ir-Mf-M^cos^ 



with the functions (|A4p given by 

V's 

V'Cb 



ilnd 



TO 



2M 
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7S 
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.-2V|'^2 



- 2Mr) sin^ 
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M - b cos ( 



(r - A^)2 - cos2 I 
1 777,2 r(r _ 2M) sin^ 9 
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i?4 



R 



1 



R = [{r - Mf + b^ ~2{r~ M)b cos 9 - M^sm^ 9]^/"^ 



(AlO) 



by virtue of the choice C — 2mA4/{b^ — Ai^) for the arbitrary costant C. It is easy to show that the linearization of 
this exact solution with respect to to is just the metric (|146|l . 



b. Superposition of two Schwarzschild black holes 



The solution corresponding to a linear superposition of two Schwarzschild black hole with masses A4 and to and 
positions z = and z = & on the z axis respectively is given by metric (jAljl with functions 

■0 = V's + -iAsb , 7 = 7s + 7Sb + 7sSb , (All) 
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where ips and 73 are given by (|A6|) . and 



7SSb 



+ Rr, 



2m 



2m 



-S(l+,2-)-E'(l-,2+) 



i?(l+_2+)-B(l-,2-) 



7Sb 



{R+ + R^f ~ im"^ 



i?9 



\/p' + (^2^)^ , 



Zj^ = z ± 7W , ^2=^ = z - (& =F to) 



(A12) 



the function 7sSb being obtained by solving Einstein's equations (jA2p . The value of arbitrary constant C can be 
determined by imposing the regularity condition (|A3I) . A unique choice of the arbitrary constant C allowing to make 
zero the function 733^ on the whole z axis does not exist: in fact it vanishes on the segment AA < z < h — m between 
the sources for C = — In ([fe^ — {M + to)^]/[6^ — {M — m)^]), and outside them (that is, for z > b + m and z < —M) 
if C is chosen to be equal to zero. 

In Boyer-Lindquist coordinates, under the transformation IjASp the metric (jAlll is given by (|A9|) with functions 



-08 = 

7sSb = 

Rf = 

= 



■ In 1 



7s 



■In 
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R^ Ri 



Ro 



R^ - 2m 
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R^ f - Am^ 
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■In 
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r - M±Mcos9 , 
{r -M)cos9±M , 



E^i±,2±) ^ r{r - 2M) sin^ d + RfR^ + Zfzf , 

R^ = [{r-Mf + (hTmf -'^{r-M){h Tm) cos 6- M^sin^ 
Z^ = {r-M) cos e-{hTm) , 



11/2 



(A13) 



by virtue of the choice C = of the arbitrary costant C . It is easy to show that also in this case the linearization of 
the solution with respect to m is just the metric ()146ll . 

APPENDIX B: THE WEYL CLASS CHARGED TWO-BODY SOLUTION 

Axisymmctric static electrovacuum solutions of the Einstein field equations can be described by the Weyl formalism 
|36|. The line element in coordinates {t,p,z,(l)) is given by Eq. ljAl|l . The electrovacuum Einstein-Maxwell field 
equations in Weyl coordinates reduce to 



2 „/,2 
P 



7,p - pI'^'^o - "0: 

7,2 - 2p[?/',p'0,z 



-2^r 



where 



X. 



X . 



and the 4-potential is determined by the electrostatic potential V only 

= -~V{p,z)dt . 



(Bl) 



(B2) 



(B3) 



Thus, once the first two equations are solved for tp and V the last two equations serve to determine 7. 

The solutions belonging to the Weyl class are characterized by the metric function ip which is a function of the 
electrostatic potential, i.e. ip = so that the gravitational and electrostatic equipotential surfaces overlap. Weyl 

[Sfij l showed that for asymptotically flat boundary conditions the unique functional relationship between ip and V is 
given by 



o2V 



Qtot 



(B4) 
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where Mtot and Qtot are the total mass and charge of the swtem, respectively. The solution can thus be written in 
terms of only one function, say f{p, z), as follows (see e.g. jSTj ) 



1 



In 



1) 



V = a 



(B5) 



where the parameter a is defined by 



(B6) 



The solution representing two Reissner-Nordstrom black holes separated by a distance h is given by [2^ 



Rt + i?r + 2Li J \Rt + R2 



2L, 



2L, 



Rt 



i?9 



(B7) 



The constants Li^2 are the lenghts of the Weyl rods generating the solution, while the parameter a is related to the 
total mass Aftot = Mi + M2 and charge Qtot = Qi + Q2 of the system through Eq. (|B6p : in addition, one has 



Ml 
M2 



1 + 0^ 

1 + a2 



Li 
L2 



so that the constants Li and L2 turn out to be given by 



Qi 
Q2 



2a 



I -a? 

2a 
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(B8) 
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(B9) 



The metric function -0 and the electrostatic potential V are found through Eq. (jB5|l . while the metric function 7 can 
be obtained by solving the field equations (|B1() 



7 = iln 



{Rt + RiY-ALl\ ( {Rt + R2r - ^Ll\ (E^i^^2-)E(i-,2 



E{1±^2± 



RfR^ 



4:R+R^ 

+ z^zt 



4i?^i?2 



^(l+,2+)-£^(l-,2-) 



c 



(BIO) 



The value of arbitrary constant C can be determined by imposing the elementary flatness condition ljA3|l We have 
that a unique choice of C allowing to make zero the function 7 on the whole z axis does not exist: in fact, it vanishes 
on the segment Li < z < b — L2 between the sources for C = — In ([6^ — {Li + L2)^]/[6^ — (Li — L2)'^]), and outside 
them (that is, for z > b + L2 and z < —Li) if C is chosen to be equal to zero. The value of 7 on the symmetry axis 
between the sources is thus given by 



lim 7 = In 



b-'-{Li+L2) 
62 - {Li - L2) 



(Bll) 



and can be made vanishing only by imposing that L1L2 = 0, which implies that both sources are critically charged, 
from Eq. l|59|) . 

The relation with standard Boyer-Lindquist coordinates in the Reissner-Nordstrom case is simply 



p = r/(r)^/^ sine* , z={r-M)cost 



(B12) 
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